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Chapter  1.  Introduction 


The  project  required  a  research  team  to  investigate  and  study  plasma  bubbles  through  their 
spread  echo  in  ionograms  limiting  to  equatorial  ionosphere,  for  the  reasons:  (1)  part  of  the 
Philippines  lies  at  the  magnetic  equator,  (2)  equatorial  ionosphere  is  not  well  understood,  (3) 
plasma  bubbles  form  in  the  equatorial  ionosphere,  and  (4)  the  equations  are  simpler  because  the 
magnetic  field  is  nearly  horizontal,  so  that  the  ionospheric  layers  are  also  nearly  horizontal.  Also 
study  plasma  bubbles  to  determine  from  an  ionogram,  the  electron  density  as  a  function  of  true 
height  which  is  not  a  very  easy  problem.  Generate  the  ionogram  as  a  plot  of  the  virtual  height  of 
the  ionosphere  as  a  function  of  the  carrier  frequency  of  the  radio  wave.  This  study  shall  generate 
an  atlas  of  all  possible  electron  density  profiles  and  derive  the  ionograms.  Using  experimentally 
obtained  ionograms  and  comparing  with  the  theoretical  ionograms,  the  corresponding  electron 
density  profile  of  the  ionogram  shall  be  found. 

This  report  consists  of  two  papers  that  have  been  submitted  for  publication.  They  constitute 
Chapters  2  (Polarization  ellipse  and  Stokes  parameters  in  geometric  algebra)  and  3  (Refractive 
index  and  energy-momentum  of  elliptically  polarized  radio  waves  in  plasma  via  geometric 
algebra). 

A  third  paper  will  be  submitted  later: 

Clint  G.  Bennett,  Quirino  M.  Sugon  Jr.,  and  Daniel  J.  McNamara,  "Simulation  of  electron 
density  profiles  and  their  ionograms  for  vertical  incidence  radio  wave  propagation  in  equatorial 
ionosphere." 
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In  this  paper,  we  use  geometric  algebra  to  describe  the  polarization  ellipse  and  Stokes  parameters.  We  show  that  a 
solution  to  Maxwell’s  equation  is  a  product  of  a  complex  basis  vector  in  Jackson  and  a  linear  combination  of  plane 
wave  functions.  We  convert  both  the  amplitudes  and  the  wave  function  arguments  from  complex  scalars  to 
complex  vectors.  This  conversion  allows  us  to  separate  the  electric  field  vector  and  the  imaginary  magnetic  field 
vector,  because  exponentials  of  imaginary  scalars  convert  vectors  to  imaginary  vectors  and  vice  versa,  while  ex¬ 
ponentials  of  imaginary  vectors  only  rotate  the  vector  or  imaginary  vector  they  are  multiplied  to.  We  convert  this 
expression  for  polarized  light  into  two  other  representations:  the  Cartesian  representation  and  the  rotated  ellipse 
representation.  We  compute  the  conversion  relations  among  the  representation  parameters  and  their  correspond¬ 
ing  Stokes  parameters.  And  finally,  we  propose  a  set  of  geometric  relations  between  the  electric  and  magnetic 
fields  that  satisfy  an  equation  similar  to  the  Poincare  sphere  equation.  ©  2011  Optical  Society  of  America 
OCIS  codes:  120.2130,  260.5430,  260.2110,  350.7420. 


1.  INTRODUCTION 

There  are  many  formalisms  to  describe  polarized  light:  Stokes 
parameters  [1,2],  Jones  calculus  [3,4],  Poincare  sphere  [5], 
Mueller  matrices  [6],  Pauli  fonnalism  [7-9],  Jackson’s  com¬ 
plex  vectors  [10],  complex  numbers  [11,12],  and  the  more 
recent  pure  operatorial  Pauli  algebraic  approach  [13-15], 
All  these  formalisms  are  dependent  on  each  other,  and  there 
are  fundamental  textbooks  [IT,  17]  in  polarization  theory  that 
establish  the  correspondences  between  them. 

Many  of  these  formalisms,  however,  use  the  imaginary 
number  simply  as  a  mathematical  artifice  for  deriving  equa¬ 
tions:  only  the  real  part  is  physically  meaningful;  the  imaginary 
part  is  discarded.  Also,  the  magnetic  field  is  not  explicitly  in¬ 
cluded  in  the  description  of  polarized  light.  In  this  paper,  we 
wish  to  show  that  imaginary  numbers  are  also  physically 
meaningful,  especially  in  the  description  of  electromagnetic 
waves  and  polarized  light.  We  shall  do  this  using  geometric 
algebra. 

Clifford  (geometric)  algebra  was  invented  by  Clifford  in  1878 
[181,  but  it  was  only  applied  to  polarized  light  by  Hestenes  in 
1971  [19],  Baylis  in  1993  [8],  and  Doran  and  Lasenby  in  2003 
[9].  Clifford  algebra  is  compact  and  coherent.  For  example, 
the  geometric  product  of  two  vectors  may  be  expressed  as  a 
simi  of  their  dot  and  cross  products  [20]: 

ab  =  a  ■  b  +  ia  x  b.  (1) 

where  i  is  a  geometric  product  of  three  orthonormal  vectors. 
Except  for  the  absence  of  the  Pauli  a  matrices,  the  geometric 
product  in  Eq.  (1)  is  similar  to  the  Pauli  identity  in  quantum 
mechanics  [21], 

Using  the  identity  in  Eq.  (1),  we  can  define  an  electromag¬ 
netic  field  £  as  a  sum  of  its  electric  field  vector  and  its 
imaginary  magnetic  field  vector: 

1084-7529/12/010089-10$15.00/0 


E  =  E  +  if  H,  (2) 

where  f  =  y/p 0/e0  is  the  ratio  of  E  to  H  for  the  radiation  field 
as  given  in  Haus  and  Melcher  [22].  This  expression  lets  us  de¬ 
fine  an  elliptically  polarized  electromagnetic  wave  in  terms  of 
the  exponential  plane  wave  functions  ei 0"t-kz): 

E  +  if  H  =  (d  +  ie2)[a+e-i('n!-te+a+)  +  a_ei(-mt~kz-s-'>].  (3) 

In  Baylis  [8],  Hestenes  [19],  and  in  Doran  and  Lasenby  [9],  the 
electromagnetic  field  is  given  F  =  E  +  iB,  because  the  mag¬ 
netic  flux  density  B  and  not  the  magnetic  field  H  is  the  phy¬ 
sical  magnetic  field  in  the  Lorentz  force.  But  we  are  following 
the  convention  of  Jancewicz  [23]  wherein  the  electromagnetic 
field  is  given  by  v/cE  +  V//7H.  Dividing  this  by  y/e  would  yield 
our  expression  in  Eq.  (2). 

We  shall  not  go  into  the  Stokes  subspace  of  the  Poincare 
sphere  [8,24],  which  is  spanned  by  the  2x2  Pauli  spin 
matrices  04  =  ej ,  rr2  =  e2,  and  <r3  =  e:l.  These  matrices,  which 
also  satisfy  the  orthonormality  relation  of  Clifford  algebra, 
act  on  an  abstract  column  vector  defined  as  the  spinor 
(a,  erlA' ,  a_e',]  ).  Note  that  despite  their  formal  similarity,  the 
Pauli  spin  matrices  tire  not  equivalent  to  the  vectors  e1;  e2,  and 
e3  that  span  the  real  physical  space. 

Doran  and  Lasenby  [9]  have  defined  polarized  light  as 

E  =  ReffCiUj  +  c2o2)e-ikx)  =  Re(£).  (4) 

We  shall  instead  use  the  authors’  exponential  relations  to 
transform  the  exponents  in  Eq.  (2)  from  i  to  -ie3.  This  allows 
us  to  conveniently  separate  the  vector  and  imaginary  vector 
parts  to  yield  the  relations  for  the  electric  and  magnetic  fields: 

E  =  e1[a+e^0,t-kz+s+'l  +  (5a) 

©  2012  Optical  Society  of  America 
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f  H  =  e2[a+eie3<w‘-fe+'5+>  +  (5b) 

In  this  form,  both  the  electric  and  magnetic  fields  are  ex¬ 
pressed  as  sums  of  counter-rotating  vectors  with  radii  a+ 
and  a_,  initial  rotation  angles  <5 .  and  <5_,  and  rotation  param¬ 
eters  t  and  z  for  time  and  space.  That  is,  the  exponential  of  an 
imaginary  vector  e±“3(i ot-kz±s±)  acfS  as  a  rotation  operator  to 
the  vectors  e,  and  e2.  Using  the  rules  of  geometric  algebra,  we 
shall  relate  this  form  with  two  other  known  representations: 
(i)  Cartesian  representation  whose  component  oscillations 
are  along  the  x  and  y  axes  and  (ii)  rotated  ellipse  representa¬ 
tion  whose  component  oscillations  are  along  the  semimajor 
and  semiminor  axes  of  the  polarization  ellipse. 

The  Stokes  parameters  S0,  Si,  S2,  and  S3  are  defined  as 
four  observable  quantities  that  characterize  the  polarization 
state  of  light  [25].  In  the  literature  [3,26],  there  is  no  adequate 
summary  of  the  Stokes  parameters  in  the  three  representa¬ 
tions  of  polarized  light,  so  we  shall  present  this  summary  here. 
We  shall  show  that  the  Stokes  parameters  are  related  to  the 
geometry  of  the  polarization  ellipse  and  that  the  first  Stokes 
parameter  S0  corresponds  to  the  magnitudes  of  the  time- 
averaged  energy  and  momentum  densities.  Furthermore,  we 
shall  define  the  field  relation  parameters: 


s0  =  E2  +  C2H2, 

(6a) 

Si  =  E2  -  C2 H2, 

(6b) 

is2  =  2ifE  •  H  =  2if|E||H|  cos  a, 

(6c) 

s3  =  2fE  x  H  =  2e3C|E||H|  sina. 

(6d) 

where  a  is  the  angle  between  the  vectors  E  and  H.  Although 
these  parameters  are  not  equivalent  to  the  classical  Stokes 
parameters,  they  are  physically  important  because  they  define 
the  relationships  between  the  electric  and  magnetic  fields  that 
are  not  found  in  the  classical  Stokes  parameters.  Moreover, 
the  field  relation  parameters  satisfy  an  equation  similar  to  that 

of  the  Poincare  sphere: 

<?2  —  <?2  1  o2  1  o2 

*0  ^  ^1  '  ^2  '  a3* 

(?) 

This  paper  is  divided  into  six  sections.  Section  1  is  the  in¬ 
troduction.  In  Section  2,  we  shall  discuss  geometric  products 
of  vectors  and  rotation  operators  in  geometric  algebra  C/:10. 
These  operators  are  expressed  in  terms  of  the  exponentials 
of  imaginary  numbers  and  of  imaginary  vectors.  In  Section  3, 
we  discuss  Maxwell’s  equation  and  its  electromagnetic  wave 
solutions  expressed  as  a  sum  of  products  of  the  circularly  po¬ 
larized  basis  and  wave  functions.  We  express  in  geometric  al¬ 
gebra  terms  three  classical  representations  of  the  polarization 
ellipse:  Cartesian  [16],  circular  [F7],  and  rotated  ellipse  repre¬ 
sentations.  In  Section  4,  we  shall  present  the  classical  Stokes 
parameters  in  all  three  representations.  We  shall  also  propose 
a  set  of  relations  between  the  electric  and  magnetic  fields  in 
terms  of  four  parameters:  the  smn  and  difference  of  their  mag¬ 
nitudes  and  their  vector  dot  and  cross  products.  We  shall 
show  that  these  magnitudes  are  related  to  the  electromagnetic 
energy  and  momentum  densities  and  the  orthogonality  and 


magnitude  relations  between  the  electric  and  magnetic  fields. 
Section  5  is  the  conclusion. 

2.  GEOMETRIC  ALGEBRA 

In  this  section  we  shall  introduce  the  products  of  basis  vectors 
and  complex  numbers  in  geometric  algebra  and  use  these  to 
derive  the  associative  geometric  product  of  vectors  as  a  sum 
of  their  dot  and  imaginary  cross  products.  By  doing  so,  we 
shall  show  that  a  general  form  of  an  element  in  geometric  al¬ 
gebra  for  three-dimensional  space  consists  of  a  scalar,  vector, 
imaginary  vector,  and  an  imaginary  scalar,  which  corresponds 
to  a  magnitude,  an  oriented  line  segment,  an  oriented  plane, 
and  an  oriented  volume,  respectively.  After  this,  we  shall  de¬ 
velop  the  formalism  for  vector  rotations.  We  shall  start  with 
exponentials  of  imaginary  scalars  as  rotation  operators  of 
complex  vectors.  Then  we  shall  proceed  with  exponentials 
of  imaginary  vectors  as  rotation  operators  of  vectors. 

A.  Vectors  and  Imaginary  Numbers 

There  are  many  vector  spaces  in  which  geometric  algebras 
can  be  generated,  but  we  shall  focus  on  the  geometric  algebra 
CI3  o,  also  known  as  the  Pauli  algebra  [23],  which  is  generated 
by  the  three  orthonormal  vectors  e1;  e2,  and  e3  of  the  physical 
space.  Any  product  of  these  three  vectors  satisfy  the  following 
orthonormality  axiom: 

£/£/,'  Cfcfy  25y/. .  (8) 

This  means  that 

e|  =  eg  =  1.  (9a) 

ejek  =  -ekej,  j  *  k.  (9b) 

In  other  words,  parallel  vectors  commute  while  perpendicular 
vectors  anticommute. 

The  algebra  also  defines  a  trivector  i, 


i  =  e!e2e3,  (10) 

which  has  the  following  properties: 

i2  =  -1,  (11) 

it  1  =  ei?.  =  e2e3,  (12a) 

ie2  =  t2i  =  e:ie, .  (12b) 

ie.3  =  e3i  =  e,e2.  (12c) 


Equation  (11)  states  that  the  trivector  i  is  an  imaginary 
number.  Equations  (12a)-(12c),  on  the  other  hand,  state  that 
i  not  only  commutes  with  vectors  but  also  transforms  these 
vectors  into  oriented  planes  perpendicular  to  them.  These 
planes  are  called  bivectors. 
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B.  Vector  Products  and  Multivectors 

Let  vectors  a  and  b  be  defined  as 

a  =  +  a2e2  +  a3e3,  (13a) 

b  =  i>iei  +  t»2e2  +  ^3e3-  (13b) 

Their  product  is 

ab  =  a  •  b  +  aAb,  ( 14) 

where 

a-b  =  a1b1+a2b2  +  a3b3,  (15a) 

aAb  =  e1e2(a1b2  -  0361)  +  e2e3(a263  -  a3b2) 

+  e3e1(a361  -aYb3).  (15b) 

Notice  that  the  product  of  the  two  vectors  a  and  b  is  the  sum  of 
a  scalar  dot  product  a  •  b  and  a  bivector  wedge  product  aAb. 
Note  that  a  •  b  is  symmetric  and  aAb  is  antisymmetric. 

To  relate  the  wedge  product  aAb  to  the  more  familiar 
cross  product  a  x  b,  we  use  the  relations  in  Eqs.  (12a)-(12c) 
to  obtain 

aAb  =  ia  x  b.  (16) 

Here,  the  bivector  (oriented  plane)  aAb  is  rewritten  in  terms 
of  the  pseudovector  a  x  b  normal  to  the  plane  formed  by  a 
and  b. 

Geometric  algebra  allows  elements  belonging  to  different 
vector  spaces  to  be  added  up;  e.g.,  the  geometric  product 
is  the  sum  of  a  scalar  and  an  imaginary  vector.  This  means 
that  the  most  general  form  of  the  sum  can  be  constructed 
by  having  a  linear  combination  of  scalars,  vectors,  imaginary 
vectors,  and  imaginary  scalars: 

A  =  a0  +  ai  +  ia2  +  ia-3,  ( 17) 

where  a0  and  a3  are  scalars  and  a1  and  a2  are  vectors.  The 
spatial  inverse  A*  of  .1  is  obtained  by  reversing  the  sign  of 
every  vector  in  A: 

A'  =  a0  -  ax  +  ia2  -  ia3,  (18) 

which  is  similar  to  the  parity  operation  in  Sakurai  [27]  except 
that  we  define  i  =  626265,  which  changes  sign  under  the  parity 
operation.  Notice  that  scalars  and  imaginary  vectors  do  not 
change  sign  under  the  parity  operation,  while  vectors  and  ima¬ 
ginary  scalars  do.  Note  that  the  dagger  notation  used  here  is 
not  the  same  as  the  reversion  operator,  such  as  in  Hestenes 
and  Sobczyk  [20],  which  reverses  the  order  of  the  vector  fac¬ 
tors,  so  that  the  reverse  of  A  is  a0  -  ai  -  ia2  -  ia3.  This  is  not 
the  same  as  Eq.  (18).  The  definition  of  spatial  inversion  in 
Eq.  (18)  is  actually  the  same  as  the  automorphic  grade  involu¬ 
tion  in  Baylis  [8]. 
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C.  Vector  Rotation  via  Exponentials  of  Trivectors 

Because  i2  =  - 1,  then  by  Euler’s  theorem, 

e±i0  _  cos  q  _j_  i  sjn  0  ( 19) 

which  is  a  sum  of  a  scalar  and  a  trivector.  From  complex 
analysis,  we  know  that 

cos  8  =  g  (e10  +  e~%e) ,  (20a) 

sin8  = -^(e18  -  e~lS).  (20b) 

To  rotate  vectors  using  the  exponential  e1’",  we  cannot  use 
ordinary  vectors  as  operands.  Instead,  we  use  the  complex 
vector  in  Jackson  [10]: 

e+  =  ei  +  ie2,  (21) 

and  we  extract  the  real  part  later.  That  is,  if  we  define 

ei+  +  ie2±  =  (er  +  it2)e±w,  (22) 

and  we  use  Eq.  (19)  in  Eq.  (22),  we  obtain 

elT  =  ex  cosd+e2  sind,  (23a) 

e2±  =  ±ei  sin  8  +  e2  cos  8,  (23b) 

after  separating  the  scalar  and  bivector  parts.  Geometrically, 
ei+  is  a  rotation  of  ei  about  e3  counterclockwise  by  an  angle  8, 
while  ex_  is  similar  except  that  the  rotation  is  clockwise.  On 
the  other  hand,  e2+  is  a  rotation  of  e2  about  e3  clockwise  by  an 
angle  8,  while  e2_  is  similar  except  that  the  rotation  is  counter¬ 
clockwise  (Fig.  T). 

D.  Vector  Rotation  via  Exponentials  of  Bivectors 

Because  (ie3)2  =  (e,  e2)2  =  -1,  then  by  Euler’s  theorem  [28], 

e±ie3e  _  cos  q  ^  ^Cg  sjn q  _  cos  0  -j-  eie2  sjj,  q  (24) 

which  is  a  sum  of  a  scalar  and  a  bivector.  From  complex  ana¬ 
lysis,  the  relations  corresponding  to  Eqs.  (20a)  and  (20b)  are 

cos  8  =  *  (eie*8  +  e~ie°e) ,  (25a) 

sind  =  -^ie3(e“39  -  e_“3,?).  (25b) 

The  product  of  e!  and  e2  with  the  bivector  ?e3  is  also  a  rotation 
of  these  two  vectors  about  e3  counterclockwise  by  an  angle 
*/2: 

erte3  =  626262  =  e2,  (26a) 

e2re3  =  e2eie2  =  -ei,  (26b) 
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where  we  used  the  relation  in  Eq.  (12c).  In  general,  if  we  want 
to  rotate  ei  and  e2  about  e3  arbitrarily  by  an  angle  6,  we  right- 
multiply  these  two  vectors  by  e  :  ,e;|W  in  Eq.  (24),  where  the  ne¬ 
gative  sign  means  clockwise  rotation  and  the  positive  sign 
means  counterclockwise  rotation: 


eig±w30  =  6l  cost?  ±  e2  sind, 

(27a) 

e2e±lea6  =  +e!  sin  6  +  e2  cos  6. 

(27b) 

Comparing  Eqs.  (27a),  (27b),  (23a),  and  (23b), 

we  get 

e1=F  =  eie™30, 

(28a) 

e2±  =  e2e+ie3e. 

(28b) 

Adding  Eqs.  (28a)  and  (28b),  and  using  Eq.  (22),  we  arrive  at 

(ej  +  it2)e±x9  =  (ei  +  ie2)e+“361. 

(29) 

Equation  (29)  is  the  exponential  conversion  relations  for 
vector  rotations. 

The  products  of  the  unit  vectors  with  the  exponential 

are 

eig±le3S  =  g+*<=3  0Cl, 

(30a) 

e2g±ie3e  _  g+*<=3«e2i 

(30b) 

e3e±les0  =  e±tC30e3. 

(30c) 

Notice  that  ei  and  e2  change  the  sign  of  the  exponent  upon 
swapping  while  e3  will  not  change  the  sign  of  the  exponent. 

3.  ELECTROMAGNETIC  THEORY 

In  this  section,  we  show  how  the  four  Maxwell’s  equations  in 
free  space  may  be  combined  in  a  single  equation  that  states 
that  the  space-time  derivative  of  the  electromagnetic  field  is 
zero,  where  the  vector  electromagnetic  field  is  a  sum  of  a  vec¬ 
tor  electric  field  and  a  bivector  magnetic  field.  We  shall  show 
that  a  solution  to  Maxwell’s  equation  is  a  circularly  polarized 
electromagnetic  wave,  and  from  this  we  construct  an  ellipti- 
cally  polarized  wave  solution.  Using  geometric  algebra,  we 
shall  transform  the  elliptically  polarized  solution  from  an  ex¬ 
pression  containing  exponentials  of  trivectors  to  another  con¬ 
taining  exponentials  of  bivectors.  After  this,  we  convert  the 
circular  basis  representation  of  elliptically  polarized  light  to 
the  Cartesian  and  rotated  ellipse  representations. 

A.  Maxwell’s  Equation 

Maxwell’s  equation  for  electromagnetic  fields  in  free  space 
can  be  expressed  in  geometric  algebra  as  a  single  equation 
[29,30]: 


where 


E  =  E  +  if  H, 

(32a) 

a  i  d  „ 

(32b) 

—  = - h  V 

dr  cdt 

are  the  electromagnetic  field  [29,31]  and  the  space-time  deri¬ 
vative  operator,  respectively.  Distributing  the  space-time  de¬ 
rivative  operator  and  applying  the  definition  of  the  geometric 
product  to  the  nabla  operator,  we  get 


V  ■  E  =  0. 

(33a) 

-^-fVxH  =  0, 

C  dt 

(33b) 

VxE  +  ^y-0, 
c  dt 

(33c) 

V  •  H  =  0, 

(33d) 

after  separating  the  scalar,  vector,  bivector,  and  trivector 
parts.  Equations  (33a)-(33d)  are  Gauss’s  law,  Ampere’s 
law,  Faraday’s  Law,  and  the  magnetic  flux  continuity  law, 
respectively. 

B.  Circularly  and  Elliptically  Polarized  Electromagnetic 
Waves 

We  assume  the  following  electromagnetic  plane  wave  solution 
to  Maxwell’s  equation: 

E±  =  E±  +  ifH±  =  e+d±i^±1. 

(34) 

where 

e+  =  e1  +  ie2, 

(35a) 

a±  =  a±elS± , 

(35b) 

1  —  g +l(mt-kz) 

(35c) 

Here,  d±  is  the  complex  amplitude  and  i// * 1  is  the  wave  func¬ 
tion.  Note  that  though  we  use  a  complex  vector  e+  that  is  pro¬ 
portional  to  Jackson’s  e+  =  e+/\/2,  we  shall  not  use  his  other 
basis  e_  =  (e  [  -  ie2)/\/ 2;  the  complex  vector  e+  is  sufficient 
for  our  purposes. 

In  order  to  understand  the  physical  meaning  of  Eq.  (34),  we 
substitute  Eqs.  (35a)-(35c)  into  our  solution, 

E±  =  E±  +  if H±  =  (ex  +  ie2)a±e:Fl(a’t_te:F5±).  (36) 

Separating  the  vector  and  bivector  parts,  we  arrive  at 
E±  =  cos((yf  -  kz+5±)  ±  e2a±  sin(a>f  -  kz+8±),  (37a) 
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Fig.  1.  Rotation  of  the  unit  vectors  e!  and  e2  counterclockwise  by 
an  angle  8  about  the  unit  vector  e3.  The  new  vectors  are  e)  =  e1ele30 
and  e2  =  e2e“ae.  The  unit  vector  e3  is  pointing  out  of  the  page. 

fH±  =  e2a±  eos(o>i  -  kz+8±)+tla±  sin  {ait  -  kz+5±).  (37b) 

Equations  (37a)  and  (37b)  describe  circularly  polarized  elec¬ 
tric  and  magnetic  fields.  If  the  subscript  of  E  and  H  is  positive, 
then  the  polarization  is  left-circularly  polarized.  If  negative, 
the  polarization  is  right-circularly  polarized. 

To  make  the  expressions  for  these  fields  more  compact,  we 
use  the  theorem  in  Eq.  (29)  for  Eq.  (36)  to  get 

E±  =  E±  +  if H±  =  (d  +  ie2 )  a±  e*®'3  (®*-te±5±) 

=  e1a±eizit*(0,t-kz±s^  ±  ie 2a±e±ie3(«,^te±«±)_  (38) 

Separating  the  vector  and  bivector  parts,  we  obtain 

E±  =  e1a±e±^(,*-b±^,  (39a) 

f  H±  =  e2a±e±ie3^mt~kz±s±) ,  (39b) 

after  factoring  out  i  in  the  second  equation.  Notice  that  the 
electric  field  E±,  the  magnetic  field  H±,  and  the  wave  vector 
k  =  A'e.  are  mutually  perpendicular,  and  that  E±  and  H±  are 
rotating  about  k  in  time  t  for  a  fixed  position  r  =  r0e3  (Fig.  2). 
If  the  sign  of  the  exponential  argument  before  ie3  is  positive, 
then  the  rotation  is  counterclockwise  in  time;  if  negative,  then 
the  rotation  is  clockwise. 

Normally  in  the  literature,  instead  of  e  =  e,  +  ie2,  only  et 
is  used,  so  that  instead  of  the  complex  monochromatic  elec¬ 
tromagnetic  plane  wave  in  Eq.  (36),  we  have 

E±  =  e1a±e+<(“*-te:F5±).  (40) 

Separating  the  vector  and  imaginary  vector  parts,  we  get 

Re(E±)  =  ei«±  cos(o>t  -  kz+5±),  (41a) 

Im(E±)  =  +-ieia±  sin(a>t  -  kz+8±).  (41b) 

By  convention,  only  the  real  part  describes  the  electric  field 
and  the  imaginary  part  is  discarded,  so  that  Eq.  (40)  only  de¬ 
scribes  a  linearly  polarized  light.  Comparing  this  result  with 
the  circularly  polarized  wave  expression  in  Eq.  (36),  we  see 
the  conciseness  and  completeness  of  geometric  algebra  that 


Fig.  2.  Rotation  of  the  electric  field  E  and  magnetic  field  H  with  zero 
initial  phase  angle.  The  initial  field  vectors  E0  and  H0  point  along  the 
e!  and  e2  directions  and  are  rotated  by  an  angle  <p  =  mt-  kz.  The  wave 
vector  k  is  pointing  out  of  the  page. 


uses  the  imaginary  number  i  to  unify  the  electric  and  magnetic 
fields  and  convert  one  to  the  other  as  the  wave  propagates 
in  space  and  time  as  shown  by  Eqs.  (37a)  and  (37b).  Thus, 
in  geometric  algebra,  a  circularly  polarized  light  is  a  more  nat¬ 
ural  expression  for  the  polarized  monochromatic  electromag¬ 
netic  plane  wave. 

Now  that  we  have  established  that  Eq.  (34)  is  a  circularly 
polarized  electromagnetic  wave,  we  now  verify  if  it  does  in¬ 
deed  satisfies  Maxwell’s  equation  in  Eq.  (31).  Substituting 
Eq.  (34)  to  Maxwell’s  equation  in  Eq.  (31),  we  get 

=  +i  ^  -  ke^j  (d  +  iz2)a±\j/±l  =  0.  (42) 

Because  co/c  =  k  by  the  wave  condition  [29]  and 

(1  -  e3)(ei  +  te2)  =  (1  _  e3)(l  +  e3)e!  =  0,  (43) 

then  the  electromagnetic  wave  Eq.  (34)  is  indeed  a  solution  to 
Maxwell’s  equation  in  Eq.  (31). 

We  may  also  rewrite  Eq.  (42)  in  a  coordinate-free  way  as 

^  =  Kc"k)(E±  +  *CH±)=°’  (44) 

where  we  used  the  relations  in  Eq.  (34)  together  with  k  =  ke3. 
Factoring  out  the  trivector  i  in  Eq.  (44)  and  separating  the 
scalar,  vector,  bivector,  and  trivector  parts,  we  get  [30] 


o 

II 

-H 

(45a) 

-E±  +fkxH±  =  0, 
c 

(45b) 

i— (TL.  -  ik  x  E,  =  0. 
c 

(45c) 

if  k  •  H±  =  0. 

(45d) 

Because  the  electric  field  E±  and  the  magnetic  field  H±  are 
both  perpendicular  to  the  wave  vector  k  as  given  in  Eqs.  (45a) 
and  (45d),  then  we  may  rewrite  the  cross  products  of  k  with 
E±  and  H±  in  terms  of  the  geometric  product: 

E±  -  i£jjyH±  =  0,  (46a) 


fH±  +  i|^E±  =  0.  (46b) 

after  factoring  out  cu/c  =  |k|.  Notice  that  these  two  equations 
are  equal.  For  example,  if  we  left  multiply  Eq.  (46b)  by  -ik/|k| 
we  get  Eq.  (46a),  because  k2  =  |kj2. 

Using  Eq.  (46b),  we  may  rewrite  the  electromagnetic  field 
E±  as 

£±=E±  +  ifH±  =  (i  +  A)e±.  (47) 
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From  this  we  see  that  i£H±  is  the  same  as  the  bivector  or 
oriented  plane  formed  by  the  unit  vector  k/|kj  and  electric 
field  E±: 

^=me±=mae±-  (48) 

C.  Representations  of  Polarization  Ellipse 

In  general,  the  solution  to  Maxwell’s  equation  with  frequency 
at  is  a  linear  superposition  of  wave  functions  \j/~l  and  i//  [291: 

E  =  b+{a+\j/ +  (49) 

which  can  be  shown  to  be  an  expression  for  elliptically  polar¬ 
ized  light,  if  we  only  work  in  the  complex  plane  by  plotting 
only  the  curve  traced  by  the  expression  a+ ijr +d_ijr~1  as  dis¬ 
cussed  in  Azzam  and  Bashara  [16]  and  Klein  and  Furtak 
[12].  We  use  Eq.  (29)  to  transform  the  general  electromagnetic 
wave  expression  in  Eq.  (49)  to  obtain 

E  =  e+(a+\j/  +  a_ijr~l),  (50) 


Fig.  3.  Circular  basis  representation  of  the  elliptically  polarized  light 
E  as  the  sum  of  two  vectors  rotating  in  opposite  directions  with  radii 
a+  and  a_  and  initial  phase  angles  <5+  and  S_.  The  semimajor  axis  is 
E1  =  a+  +a_,  the  semiminor  axis  is  E2  =  a+  -  a_,  and  the  tilt  angle  is 

P/2  =  (5++8_)/2. 

phasors  in  complex  analysis  as  an  accessory  to  vector  analy¬ 
sis,  we  do  not  take  the  real  vector  part  and  discard  the  ima¬ 
ginary  vector  part  of  the  electric  field  E  because  the  result  of 
the  rotations  is  still  a  vector. 

Using  Euler’s  theorem  in  Eq.  (24)  and  the  identity  eri^  =  e2 
in  Eq.  (26a),  Eq.  (53)  becomes 


where 

d±=a±e“3S±,  (51a) 

^±1  =  e±^(cot-kz)  _  (51b) 

Note  that  ijr  =  in  Eq.  (49)  is  an  exponential  of  an  ima¬ 
ginary  number  or  trivector,  while  i jr  =  hr  Eq.  (50)  is 

an  exponential  of  an  imaginary  vector  or  bivector. 

Separating  the  vector  and  imaginary  vector  parts  of  the 
transformed  electromagnetic  wave  expression  in  Eq.  (50) 
yields  the  electric  and  magnetic  fields: 

E  =  ei(d+vi  +  (52a) 

CH  =  e2(a+vl  +  «_iy_1).  (52b) 


E  =  ei(a+  cos(<p  +  <5+)  +  a_  cos(</>  -  <5_)) 


+  e2(o+  sin(i/i  +  <5+)  -  a_  sin((p  -  8_)). 

(55) 

This  may  be  rewritten  as 

E  =  Ex  cos((/>  +  +  Ey  cos (</>  +  52/)e2, 

(56) 

where 

Ex  =  [aa  +  a?  +  4a+a_  cos(5+  +  <5_)]^, 

(57a) 

Ey  =  [a\  +  a?  -  4a+a_  cos(5+  +  5_)]^, 

(57b) 

(  a,  sin 5+  -  a  sin 8  \ 

8X  =  arctan  - z -  , 

\«+cos5+  +  a  cos  8  J 

(57c) 

Because  the  electric  and  magnetic  fields  have  the  same  ellip¬ 
tically  polarized  wave  expression  except  that  the  unit  vectors 
rotated  for  the  electric  field  is  e,  and  that  for  the  magnetic 
field  is  e2,  then  it  is  sufficient,  from  this  point  onward,  to  dis¬ 
cuss  only  the  electric  field. 

Using  the  definitions  of  d±  and )// 1 1  in  Eqs.  (51a)  and  (51b), 
we  may  rewrite  the  electric  field  in  Eq.  (52a)  as 

E  =  ex  [a+e“3^+a+)  +  a.e-^-S- )]  _  (53) 

where 

4>  =  cot-  kz.  (54) 


Sy  =  arctan 


a_  cos 5_  -  a,  cos 8 


I  °- 


\  a 1  sin  <5+  -I-  a_  sin  S_ 


(57d) 


Equation  (56)  is  the  Cartesian  basis  representation  of  ellipti¬ 
cally  polarized  light,  interpreted  as  a  sum  of  two  orthogonal 
oscillations  with  Ex  and  Ey  as  the  amplitudes  along  the  x  and 
y  axes,  respectively,  with  8X  and  8y  as  the  initial  phase  angles 
of  oscillation  (Fig.  4).  Equations  (64a)  and  (57d)  are  the  ex¬ 
pressions  of  the  Cartesian  basis  parameters  in  terms  of  the 
circular  basis  parameters.  Their  inverse  relations  Eire 


a±  =  l  [El+2ExEy  sin(±(5x  -  5y))  +  E$,  (58a) 


Equation  (53)  is  the  circular  basis  representation  of  elliptically 
polarized  light,  interpreted  as  a  sum  of  two  vectors  of  lengths 
a+  and  a_  that  are  rotating  counterclockwise  and  clockwise, 
respectively,  with  initial  phase  angles  <5+  and  8_  from  the  po¬ 
sitive  x  Eixis  (Fig.  3).  Notice  that  the  exponentials  e163^' 1  '5l  and 
e-ie:,i7/)-rj  1  acf  as  vector  rotation  operators  on  the  vectors  a.ei 
Eind  a._ei.  Thus,  unlike  in  Klein  and  Furtak  [12],  which  uses 


5+ 


=  arctan 


'Ex  cos  5X  -  Ey  cos  5y  \ 
Ex sin Sx  +  Ey  sin Sy  ’  J’ 


(58b) 


5_ 


-  arctan 


'  Ex  sin  5X  -  Ey  cos  Sy  \ 
Ex  cos  Sx  +  Ey  sin  Sy) 


(58c) 
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Fig.  4.  Cartesian  basis  representation  of  elliptically  polarized  light 
E5  s  as  a  sum  of  two  orthogonal  oscillations  with  amplitudes  Ex 
and  Ey,  and  initial  phase  angles  8X  and  5y. 


Fig.  5.  Rotated  ellipse  representation  of  the  elliptically  polarized 
light  Es„  as  a  rotated  ellipse  with  semimajor  axis  £),  semiminor  axis 
E2  and  tilt  angle  /?/ 2. 


On  the  other  hand,  if  we  define 

a±  =  l-(El±E2),  (59a) 

<5±  =  ±<5  +  ^ ,  (59b) 


and  substitute  these  relations  into  Eq.  (53),  we  get 


(El  +  2 ExEy  sin(4 


5y)  +  E^)  2 


+  (El  -  2 ExEy  sin(4  -  8y )  +  £$ 


(63a) 


E2 


2  (El  +  ZExEy  sin)^  -  5V )  +  )* 

-  (El  -  2 ExEy  sinfe  -  8y)  +  . 


(63b) 


E  =  ej 


*  (£)  +  E2)eie^+V  +  *  (Ex  -  E2)e~ie^ *+*) 


ei*zP! 2, 
(60) 


Using  Euler’s  theorem  in  Eq.  (24)  and  the  identity  ei?e3  =  e2  in 
Eqs.  (26a)  and  (60)  becomes 

E  =  [Ey  cos(<j>  +  <5)e!  +  E2  sin(^>  +  8)t2]el‘3^2.  (61) 

Equation  (61)  is  the  rotated  ellipse  representation  of  ellipti¬ 
cally  polarized  light,  interpreted  as  an  ellipse  with  initial  phase 
angle  8,  semimajor  axes  lengths  El  and  E2,  that  is  rotated 
counterclockwise  about  e3  by  an  angle  p/2  (Fig.  5).  The  ad¬ 
vantage  of  using  geometric  algebra  in  this  representation  is 
that  the  ellipse  El  cos(rp  +  5)ej  +  E2  sin(<j6  +  <5)e2  can  be  easily 
rotated  via  the  action  of  the  exponential  of  the  bivector  ie3. 
Equations  (59a)  and  (59b)  are  the  circular  basis  parameters 
in  terms  of  the  rotated  ellipse  representation  parameters. 
The  inverse  relations  are 


Ei  =  a+  +  a_. 

(62a) 

s' 

1 

+ 

II 

CM 

kq 

(62b) 

S  =  ^(8+-8_), 

(62c) 

(62d) 

For  completeness,  we  express  the  rotated  ellipse 
parameters  (Ei,  E2,  8,  p/2)  in  terms  of  the  Cartesian  basis 
parameters  (Ex,  Ey,  8X,  8V ): 
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1 

2 


arctan 


+  arctan 


( Ex  cos  8X  -  Ey  cos  8y\ 
v  Ex  sin  8X  +  Ey  sin  8y  J 
'  Ex  sin  8X  -  Ey  cos  8y  \ " 
Excos8x  +  Eysm8y) 


(63c) 


P 

2 


1 

2 


arctan 


-  arctan 


/ Ex  cos  8X  -  Ey  cos  8y\ 
\  Ex  sin  ST  +  Ey  sin  8y  ) 

'  Ex  sin  Sx  -  Ey  cos  8y  \ 
Ex  cos  8X  +  Ey  sin  8y)  _ 


Their  inverse  relations  are 

„  \E\  +  E\  +  (E\  -  E\)  cos/?]  I 

Ex  =  j- - 2 - J  ’ 

F  _  \E\  +  -Ef  ~  (Ef  -  -Ef)  cos/?]  i 


(63d) 


(64a) 

(64b) 


/£)  sin 8  cos/J/2  +  E2  cos <5  sin/J/2 
x  arctan  cos  ^  cog  pj2-  e2  sin  8  sin/J/2 


8y  =  arctan 


/  El  sin  8  sin  p/2  -  E2  cos<5cos/?/2 


\Ei  cos  8  sin  p/2  +  E2  sin  <5  cos  p/2 


(64c) 

(64d) 


4.  STOKES  AND  FIELD  RELATION 
PARAMETERS 

In  this  section,  we  derive  the  Stokes  parameters  in  the  rotated 
ellipse  and  circular  basis  representations  starting  with  the 
Stokes  parameters  in  the  Cartesian  basis  representation,  by 
using  the  basis  parameter  relations  that  are  given  in  the  pre¬ 
vious  section.  After  this  we  discuss  the  field  relation  param¬ 
eters  that  also  satisfy  an  equation  similar  to  the  Poincare 
sphere  equation. 
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A.  Classical  Stokes  Parameters 


The  Stokes  parameters  in  the  Cartesian  basis  representation 
in  Eq.  (56)  are  widely  known  in  the  literature.  They  are  de- 

fined  as  [11,32,33] 

S0=E*+E%, 

(65a) 

S\  =El-El, 

(65b) 

S2  =2ExEycos{Sx-Sy), 

(65c) 

S3  =  2ExEysm{8x-5y). 

(65d) 

The  parameter  S0  describes  the  total  intensity  of  the  optical 
field,  S i  describes  the  preponderance  of  linearly  horizontally 
polarized  light  over  linearly  vertically  polarized  light,  S2  de¬ 
scribes  the  preponderance  of  linear  +45°  polarized  light  over 
linear  -45°  polarized  light,  and  the  fourth  parameter  S3  de¬ 
scribes  the  preponderance  of  right-circularly  polarized  light 
over  left-circularly  polarized  light  [25].  The  definition  of  the 
Stokes  parameters  in  Eqs.  (65a)-(65d)  will  be  used  in  deriving 
the  Stokes  parameters  for  the  rotated  ellipse  and  the  circular 
basis  representations. 

We  can  show  that  the  Stokes  parameters  in  Eqs.  (65a)- 
(65d)  can  be  expressed  in  the  rotated  ellipse  representation 
(Ei,  E2,  S,  p/2)  as 


S0  =  Ef  +  El 

(66a) 

S)  =  (Ef  -El)  cos p. 

(66b) 

S2  =  (El~El)smp. 

(66c) 

S3  =  2E1E2, 

(66d) 

where  we  used  the  relations  in  Eqs.  (64a)-(64d),  together  with 
the  identities 

sini  arctan  0)  =  — -==  . . 

Vi  +  e2 

(67a) 

,  N  l 

cos(arctan^)  =  — 

Vi  +  e2 

(67b) 

Notice  that  the  Stokes  parameters  Eqs.  (66a)-(66d)  do  not  de¬ 
pend  on  the  elliptical  phase  angle  S,  but  only  on  the  geometric 
parameters  of  the  ellipse:  semimajor  axis  Elt  semiminor  axis 
E2,  and  tilt  angle  p/2  (Fig.  5),  assuming  that  Ei  is  greater 
than  E2. 

The  Stokes  parameters  in  the  circular  basis  representation 
[17,28]  can  be  derived  by  substituting  relations  Eqs.  (62a), 
(62b),  and  (62d)  in  the  Stokes  parameters  of  the  rotated 
ellipse  representation  in  Eqs.  (66a)-(66d): 


S0  =  2  (a2_  +  al), 

(68a) 

S i  =  4a+a_  cos(5+  -1-  <5_), 

(68b) 

S2  =  4 a+a_  sin(5+  +  5_), 

(68c) 

S3  =  2(a2  -  a%). 

(68d) 

Equations  (68a)-(68d)  say  that  the  Stokes  parameters  only 
depend  on  the  radii  a_  and  a+  of  the  counter-rotating  vectors 
in  the  circular  basis  and  the  sum  S+  +  S_  of  the  initial  phase 
angles. 

B.  Field  Relation  Parameters 

To  determine  the  relationship  between  E  and  H  as  defined 
in  Eqs.  (52a)  and  (52b),  we  first  take  their  squares  and  their 
products: 

E2  =  £2H2  =  a\  +  a?  +  2a+a_  cos(2</>  +  <5+  +  <5_),  (69a) 

£EH  =  -fHE  =  ie3[a2  +  a?  +  2a.+a_  cos(2<j6  +  <5+  4-  <5_)], 

(69b) 

where  we  used  the  relations  (30a)  and  (30b)  and  the  exponen¬ 
tial  definition  of  the  cosine  in  Eq.  (20a).  Adding  and  subtract¬ 
ing  these  squares  and  products,  we  get 

s0  =  E2  +  f2 H2  =  2[a2  +  al  +  2a+a._  cos(2^  +  <5+  4-  <5_)], 


(70a) 

sx  =  E2  -  C2H2  =  0, 

(70b) 

is2  =  2iCE  ■  H  =  0, 

(70c) 

s3  =  2fE  x  H  =  2 e3[a|  +  a2  +  2a+a_  cos(2 (p  +  <5+  +  <5_)] , 

(70d) 

where  s0,  Sj,  s2,  and  s3  are  what  we  shall  define  as  the  field 
relation  parameters.  Equations  (70a)  and  (70d)  are  new  re¬ 
sults  for  elliptically  polarized  light,  but  Eqs.  (70b)  and  (70c) 
are  familiar  relations  for  electromagnetic  waves.  The  first  re¬ 
lation  shows  that  the  energy  density  U  =  (e0E2  +  /i0H2)/2  is 
varying  in  space  and  time  according  to  the  wave  function  ar¬ 
gument  2 cf)  =  2 (cot  -  kz )  [8,29].  The  second  relation  shows  that 
the  magnitudes  of  the  electric  and  magnetic  fields  are  propor¬ 
tional,  |E|  =  f|H  [34],  The  third  relation  shows  that  the  dot 
product  of  the  electric  and  magnetic  fields  is  zero,  which 
means  that  E  and  H  are  perpendicular.  And  the  fourth  relation 
shows  that  the  momentum  density  S/c2  =  (E  x  H)/c2  is 
proportional  to  the  energy  density  U  =  |S|/c,  which  corre¬ 
sponds  to  the  well-known  relation  for  light  E  =  pc  [35], 
Notice  that  if  we  set  <f>  =  0,  which  means  t  =  0  at  z  =  0,  then 
the  field  relations  (70a)  and  (70d)  reduce  to 
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SoU=o  =  E2+f2H2=S0+S1.  (71a) 

S3|^o  =  2fExH  =  e3(So+Si).  (71b) 

where  S0  and  S)  are  the  classical  Stokes  parameters. 

Taking  the  time  average  over  one  optical  period  r  =  '2n/u> 
of  Eqs.  (70a)-(70d)  as  similarly  done  in  Baylis  [8],  we  get 

{s0)  =  <E2  +  f2H2)  =  2  (a2  +  a2.),  (72a) 

(sj)  =  <E2  -  f2 H2)  =  0,  (72b) 

(is2)  =  (2ifE  •  H)  =  0,  (72c) 

(s3>  =  (2fE  x  H)  =  2e3(a2  +  a2).  (72d) 


Notice  that  the  magnitudes  of  the  energy  and  momentum 
densities  are  equal  to  the  first  Stokes  parameter  S0. 

Taking  the  spatial  inverse  of  the  electromagnetic  field  E  in 
Eq.  (50),  we  obtain 

i,t  =  -  E  +  iC H,  (73) 

where  E  and  H  are  defined  in  Eqs.  (52a)  and  (52b),  respec¬ 
tively.  It  can  be  shown  that  the  field  relation  parameters 
defined  in  Eqs.  (70a)-(70d)  are  related  to  the  product  of  the 
electromagnetic  field  and  its  spatial  inverse: 


EE  =  Sl  +  is2,  (74a) 

E*&  =  -  is2,  (74b) 

&E=s  o-ss,  (74c) 

EE1  =  -s0  -  s3.  (74d) 


Notice  that  using  the  expansions  in  Eqs.  (74a)  and  (74b),  we 
may  combine  the  relations  in  Eqs.  (72b)  and  (72c)  into 
EE  =  E^E^  =  0. 

Now,  the  product  of  E‘  E  with  its  spatial  inverse  can  be 
expressed  as 

{E  Ef )  (E^E)  =  (F# )  (EE) ,  (75) 

because  the  products  are  associative  and  E'E^  is  a  scalar  by 
Eq.  (74b).  Substituting  Eqs.  (74a)-(74d)  into  Eq.  (75),  we  get 
the  relation  [31] 

s2-s2  =  s2  +  sl.  (76) 

Equation  (76)  is  similar  in  form  to  the  Poincare  sphere 
equation  [2,3]. 

Though  the  field  relation  parameters  s0,  s1;  is2,  and  s3  in 
Eqs.  (70a)-(70d)  satisfy  an  equation  similar  to  the  Poincare 
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sphere,  these  parameters  are  not  identical  to  the  Stokes 
parameters.  But  even  then,  the  field  relation  parameters  com¬ 
plement  the  Stokes  parameters  because  the  former  provides  a 
geometric  interpretation  for  the  geometric  relations  of  the 
electric  field  E  and  the  magnetic  field  H.  For  monochromatic 
plane  waves  that  we  considered  in  this  paper,  the  E  and  H  are 
orthogonal  in  direction  and  proportional  in  magnitude  to  each 
other  so  that  .S|  and  s2  are  zero.  But  when  light  is  a  mixture  of 
several  polarization  states  propagating  in  different  directions, 
the  net  electric  and  magnetic  fields  of  the  superposition  are 
not  necessarily  perpendicular  and  proportional  so  that  S] 
and  s2  will  not  be  zero.  From  these  two  parameters  we  can 
determine  the  relative  magnitudes  and  directions  of  the  time- 
averaged  electric  and  magnetic  fields. 

5.  CONCLUSION 

In  this  paper,  we  have  presented  a  geometric  algebra  formal¬ 
ism  for  polarized  light.  We  started  with  Maxwell’s  equation  in 
free  space  and  proposed  an  electromagnetic  plane  wave  solu¬ 
tion  in  terms  of  the  product  of  a  circular  basis  vector  with  the 
linear  combination  of  a  wave  function  and  its  complex  con¬ 
jugate.  We  transformed  the  complex-valued  wave  functions 
and  their  coefficients  into  terms  involving  the  exponentials 
of  an  imaginary  vector.  This  transformation  allows  us  to 
immediately  separate  the  electric  and  magnetic  fields  of  the 
electromagnetic  wave  by  taking  its  vector  and  imaginary 
vectors  parts. 

Because  the  electric  and  magnetic  fields  are  similar  in  form, 
we  only  considered  the  electric  field  component  as  done  in  the 
literature.  This  electric  field  is  the  sum  of  two  counter-rotating 
vectors  with  the  exponentials  of  imaginary  vectors  acting  as 
vector  rotation  operators.  Using  the  identities  in  geometric  al¬ 
gebra,  we  showed  that  this  representation  is  equivalent  to  the 
two  other  representations:  (i)  the  rotated  ellipse  representa¬ 
tion  whose  component  oscillations  are  along  the  semimajor 
and  semiminor  axes  of  the  polarization  ellipse  and  (ii)  Carte¬ 
sian  basis  representations  whose  component  oscillations  are 
along  the  x  and  y  axes.  We  computed  the  relations  between 
the  three  sets  of  representation  parameters.  We  also  com¬ 
puted  the  Stokes  parameters  in  the  different  representations 
and  provided  interpretations  based  on  the  geometry  of  the  el¬ 
lipse.  We  showed  that  the  first  Stokes  parameter  corresponds 
to  the  magnitudes  of  the  time-averaged  energy  and  momen¬ 
tum  densities.  We  showed  that  the  electric  and  magnetic  fields 
are  perpendicular  in  direction  and  proportional  in  magnitude. 
We  also  defined  the  field  relation  parameters  in  terms  of  the 
sum  and  difference  of  the  squares  of  the  electric  and  magnetic 
fields  and  their  dot  and  cross  products.  And  like  the  Stokes 
parameters,  they  also  satisfy  an  equation  similar  to  the 
Poincare  sphere  equation. 

The  presentation  in  this  paper  of  polarized  fight  via  geo¬ 
metric  algebra  has  several  advantages  over  the  other  formal¬ 
isms.  The  first  advantage  is  the  unifying  power  of  geometric 
algebra.  Geometric  algebra  unifies  complex  and  vector  ana¬ 
lyses,  the  four  Maxwell’s  equations,  and  the  electric  and  mag¬ 
netic  wave  solutions.  The  second  advantage  is  that  we  do  not 
disregard  imaginary  quantities  at  the  end  of  each  computation 
as  is  usually  done  in  the  literature.  Instead,  we  use  imaginary 
vectors  to  express  rotation  operators  and  to  distinguish  the 
magnetic  from  the  electric  field.  The  third  advantage  is  we 
can  use  the  machinery  of  geometric  algebra  to  convert 
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different  representations  of  polarized  light.  The  exponential 
rotation  operators  of  vectors  are  used  to  describe  the  tilt  of 
the  polarization  ellipse  and  to  describe  its  counter-rotating 
vector  components.  The  fourth  advantage  is  that  we  define 
the  field  relation  parameters  which  describe  the  geometric  re¬ 
lations  between  the  electric  and  magnetic  fields.  The  vanish¬ 
ing  of  the  second  and  third  field  relation  parameters  means 
that  the  electric  and  magnetic  fields  are  perpendicular  to 
each  other  and  proportional  in  magnitude.  The  nonvanishing 
first  and  fourth  field  relation  parameters  correspond  to  the 
electromagnetic  energy  and  momentum  densities,  whose  time 
averages  are  proportional  to  the  first  Stokes  parameter. 

In  a  future  work,  we  shall  extend  our  geometric  algebra 
formalism  to  optical  elements  such  as  rotators,  attenuators, 
phase  shifters,  polarizers,  and  birefringent  crystals,  and  dis¬ 
cuss  the  significance  of  the  field  relation  parameters.  We  shall 
also  discuss  partially  polarized  light  and  its  relationships  to 
Fourier  analysis  and  Hilbert  spaces  in  geometric  algebra. 
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Abstract 

In  this  paper,  we  revisit  the  problem  of  imaginary 
refractive  index  in  unmagnetized  plasma  for  ellipti¬ 
cally  polarized  light.  We  verified  Budden’s  statement 
that  the  electric  and  magnetic  field  of  the  electromag¬ 
netic  wave  in  plasma  would  be  parallel,  so  that  the 
Poynting  vector  is  zero.  But  we  also  showed  that 
the  energy-momentum  of  the  wave  is  not  zero:  the 
energy  density  of  the  wave  is  negative  and  its  mo¬ 
mentum  density  is  opposite  to  that  if  the  wave  were 
in  a  medium  with  a  real  refractive  index.  The  for¬ 
malism  that  use  is  the  Clifford  (geometric)  algebra 
C/3,0,  also  known  as  the  Pauli  algebra. 

1  Introduction 

It  is  known  that  for  an  isotropic,  unmagnetized  iono¬ 
sphere,  the  square  of  the  refractive  index  n  for  an 
electromagnetic  wave  with  angular  frequency  u>  is 
given  by  [1] 


where  N  is  the  ionosphere’s  electron  density,  q  is  the 
electron’s  charge,  m  is  the  electron’s  mass,  and  eo  is 
the  permittivity  of  free  space. 

By  the  law  of  Trichotomy,  there  are  only  three  pos¬ 
sibilities:  n 2  >  0,  n2  =  0,  and  n 2  <  0.  The  in¬ 
terpretation  of  the  first  two  cases  is  straightforward, 
but  not  the  third.  Budden  analyzed  the  relations 
nEy  =  — TLX  and  nEx  =  'Ey  and  concluded  that  these 
hold  whether  n  is  real  or  complex)] Budden  1966  pp 
38  41.  Furthermore,  according  to  him: 


If  n2  is  negative  and  real,  then  n  is  purely 
imaginary....  This  appears  to  represent  a 
wave  travelling  with  infinite  velocity.  Ev¬ 
ery  Held  component  varies  harmonically  in 
time,  but  there  is  no  harmonic  variation  in 
space.  The  phase  is  the  same  for  all  values 
of  z.  A  disturbance  of  this  kind  is  called  an 
‘evanescent’  wave....  In  an  evanescent  wave 
n  is  purely  imaginary,  so  that  the  electric 
and  magnetic  fields  are  in  quadrature.  This 
means  that  the  time  average  value  of  the 
Poynting  vector  n  =  E  x  H  is  zero.  Hence, 
there  is  no  net  energy-flow  for  such  a  wave 
in  an  isotropic  medium.  For  an  evanescent 
wave  in  an  evanescent  medium,  there  can 
be  some  flow  of  energy  perpendicular  to  the 
wave- normal.  [3] 

In  this  paper,  we  shall  revisit  the  problem  of  en¬ 
ergy  flow  for  imaginary  refractive  index.  We  shall 
show  that  even  if  its  the  Poynting  vector  is  zero,  an 
elliptically  polarized  electromagnetic  wave  would  still 
possess  energy  and  momentum,  except  that  the  en¬ 
ergy  density  would  be  negative  and  the  momentum 
density  would  be  opposite  to  that  if  the  refractive  in¬ 
dex  were  a  real  number.  We  shall  also  show  that  for 
n  =  0,  the  energy  and  momentum  densities  of  the 
wave  would  be  zero. 

We  shall  divide  this  paper  into  six  sections.  Sec¬ 
tion  1  is  Introduction.  In  Section  2,  we  shall  discuss 
axioms  and  theorems  of  Clifford  (geometric)  algebra 
C/3,0,  such  as  the  orthonormality  axiom,  the  Pauli 
identity,  and  spatial  inversion.  We  shall  also  discuss 
exponential  functions  and  their  applications  in  vector 
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rotations  in  the  plane.  In  Section  3,  we  shall  sum¬ 
marize  Classical  Electrodynamics  by  discussing  the 
Maxwell’s  equation,  the  electromagnetic  wave  equa¬ 
tion,  and  the  electromagnetic  energy-momentum  .  In 
Section  4,  we  shall  derive  the  electromagnetic  wave 
equation  in  plasma  by  computing  the  currents  due  to 
an  elliptically  polarized  electric  field.  We  then  deter¬ 
mine  the  electromagnetic  wave  equation  for  plasma 
and  compute  the  refractive  index.  We  shall  find  elec¬ 
tromagnetic  wave  solutions  to  this  plasma  wave  equa¬ 
tion  and  compute  their  energy-momentum  for  differ¬ 
ent  cases  of  refractive  indices.  Section  5  is  Conclu¬ 
sions. 

2  Geometric  Algebra 

2.1  Vectors  and  Imaginary  Numbers 

The  Clifford  (Geometric)  algebra  Cl3fi  is  generated 
by  three  vectors  ei,  e2,  and  e3  which  satisfy  the  or¬ 
thonormality  axiom  [4] 

&j^k  T  ^ jk  ■  (2) 

In  other  words, 

e|  =  1,  (3a) 

=  j  k'  (3b) 

In  other  words,  the  square  of  a  unit  vector  is  unity, 
while  the  product  of  two  perpendicular  vectors  anti¬ 
commute. 

Let  us  define  the  unit  trivector  i : 

i  =  e  ie2e3.  (4) 

We  can  show  that  the  trivector  i  has  the  following 
properties: 


z2  =  - 

■1 

(5) 

ie  i 

=  ^2^3 

=  eii, 

(6a) 

*e2 

=  e3ei 

=  e2i, 

(6b) 

*e3 

=  ele2 

=  e3i, 

(6c) 

(6d) 

That  is,  the  trivector  i  is  an  imaginary  number  that 
commutes  with  all  vectors. 

Let  a  and  b  be  two  vectors  spanned  by  ei,  e2,  and 
e3.  We  can  show  that  the  geometric  product  of  a  and 
b  is  given  by  the  Pauli  identity  [5] 

ab  =  a  •  b  +  i( a  x  b),  (7) 

where  ab  and  ax  b  are  the  dot  and  cross  products  of 
the  two  vectors.  From  the  definition  of  the  geometric 
product,  it  is  easy  to  see  that 


a  •  b  =  -(ab  +  ba), 

(8a) 

a  x  b  =  —  (ab  —  ba). 

(8b) 

Thus,  the  dot  product  of  two  vectors  is  half  their 
anticommutator,  while  their  imaginary  cross  product 
is  half  their  commutator. 

In  general,  any  object  A  in  Cl3to  may  be  expressed 
as 

A  =  Aq  +  Ai  +  iA.2  +  2A3,  (9) 

which  we  refer  to  as  a  cliffor.  The  spatial  inverse 
AA  is  obtained  by  replacing  each  vector  in  AA  by  its 
negative: 

AA  =  A0  —  Ai  +  iA2  —  iA3.  (10) 

Notice  that  under  spatial  inversion,  only  the  vector 
and  trivector  parts  change  sign;  the  scalar  and  bivec¬ 
tor  parts  remain  unchanged. 

Let  A  and  B  be  two  cliffors  in  Cl3fi.  The  spatial 
inverse  of  the  sum  and  product  of  these  two  cliffors 
are  given  by  [6] 

(A  +  B)*  =A*  +&,  (11a) 

(AB)*  =  A*B*.  (lib) 

In  other  words,  the  spatial  inversion  operator  (  *)  is 
distributive  over  addition  and  multiplication.  De¬ 
spite  their  similarity  in  appearance,  the  spatial  in¬ 
version  operator  is  not  the  same  as  the  reversion  op¬ 
erator  used  in  the  literature,  which  interchanges  the 
order  of  vector  factors  in  a  cliffor:  the  reversion  op¬ 
erator  changes  the  sign  of  the  bivector  and  trivec¬ 
tor  parts,  but  leaves  the  scalar  and  vector  parts  un¬ 
changed. 


2 


14 


2.2  Vector  Rotations  and  Exponential 


Functions 

Since  (ie3)2  =  —  1,  then 

ele3d  =  cos  9  +  ze3  sin  9.  (12) 

Left-  and  right-multiplying  Eq.  (12)  with  ei,  e2,  or 
e.3,  we  obtain 

eiele:}ff  =  ei  cos  9  +  e2  sin#  =  e~tB30ei,  (13a) 
e2ele361  =  — ei  sin#  +  e2  cos#  =  e_le3Se2,  (13b) 
e3ele30  =  e3  cos  #  +  sin#  =  e'e3l9e3.  (13c) 


Geometrically,  eie*®3®  and  e2ele3®  are  rotations  of  e3 
and  e2  about  e3  counterclockwise  by  an  angle  #. 

Let  us  define  the  vector  s  on  the  xy—  plane  as 

s  =  eisi  +  e2s2  =  ei  se1*30 .  (14) 


The  product  of  s  and  s'  is 


(20) 


Using  the  Pauli  identity  in  Eq.  (7)  together  with  the 
Euler’s  formula  in  Eq.  (12),  and  separating  the  scalar 
and  bivector  parts  of  Eq.  (20),  we  arrive  at 

s  •  s'  =  si  si  +  s2s(,  =  ss'  cos( — 9  +  #'),  (21a) 

s  x  s'  =  ie3(s!si  —  s'1s2)  =  ie3ss'  sin(— 9  +  9'), 

(21b) 


which  are  the  familiar  relations  for  the  dot  and  cross 
products  of  two  dimensional  vectors  in  rectangular 
and  polar  forms. 


3  Classical  Electrodynamics 


From  this  we  obtain 


3.1  Maxwell’s  Equation 


In  a  medium  characterized  by  dielectric  permittiv¬ 
ity  eo  and  magnetic  permeability  yo,  the  Maxwell’s 
9  =  tan_1(s2/si),  (15b)  equation  in  electrodynamics  is  given  by 


and 


si  =  s  cos  6,  (16a) 

s2  =  ssin#.  (16b) 

Equations  (15a)  to  (16b)  are  the  familiar  transfor¬ 
mation  relations  for  the  polar  and  rectangular  coor¬ 

dinates. 

The  vector  s  in  Eq.  (14)  may  also  be  written  as[7] 
s  =  ds  =  s*e1;  (17) 

where 

s  =  si  +  ie3s2  =  seie30 ,  (18a) 

s*  =  si  —  *e3s2  =  se~ie30.  (18b) 

We  say  that  s  is  the  complex  representation  of  the 
vector  s  and  s*  is  the  complex  conjugate  of  s. 

Let  s'  be  another  vector  in  the  xy— plane: 

s  =  eisi  +  e2s2  =  eise'®3^.  (19) 


dE 

dr 

where 

d_  _  1 
dr  c 
E  =  E 

J  =  P' 


Co/,  (22) 


d  r, 

at  +  v’ 

(23a) 

+  *CoH, 

(23b) 

:  +  j- 

(23c) 

Note  that  the  speed  of  light  c  =  l/^/yoCo  and  the 
intrinsic  impedance  £0  =  yVo  Ao  ■  In  words,  the 
Maxwell’s  equation  states  that  the  spacetime  deriva¬ 
tive  d/df  of  the  electromagnetic  field  E  is  propor¬ 
tional  to  the  spatial  inverse  of  the  charge-current  den¬ 
sity  3. 

Let  us  expand  the  Maxwell’s  Equation  in 
Eq.  (22):  [8,  9] 

(i|+v)(E  +  iCoH)=Co(pc-j).  (24) 
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Using  the  Pauli  identity  in  Eq.  (7)  and  separating  the 
scalar,  vector,  bivector,  and  trivector  parts,  we  arrive 
at 


for  the  electric  and  magnetic  field  intensities  [11],  re¬ 
spectively. 

If  the  charge  distribution  p  is  nearly  constant,  then 


V  •  E  =  — ,  (25a) 

e0 

1  BT*' 

-CoVxH  =  -Coj,  (25b) 

c  at 

au 

Mo  ~rr  +  V  x  E  =  0,  (25c) 

at 

CoV  •  H  =  0,  (25d) 

which  are  the  Gauss’s  law,  Ampere’s  law,  Faraday’s 
law,  and  magnetic  flux  continuity  law,  respectively. 


Vp  =  0,  (31) 

so  that  the  wave  equation  in  Eq.  (30b)  simplifies  to 

□E  =  -Mo§-  (32) 

Later  we  shall  express  the  current  density  j  in  terms 
of  the  electric  field  intensity  E,  in  order  to  deter¬ 
mine  the  refractive  index  and  energy-momentum  of 
the  electromagnetic  wave  in  plasma. 


3.2  Wave  Equation 


The  spatial  inverse  of  the  spacetime  derivative  oper¬ 
ator  is 


d  Id 
i  c  dt 


(26) 


Applying  this  operator  to  the  Maxwell’s  equation  in 
Eq.  (22),  we  obtain 


OE  = 


where 


□  = 


d2 

df^df 


1  d2 
c 2  dt2 


is  the  d’Alembertian  operator. 
Expanding  Eq.  (27), 


V2 


(27) 


(28) 


□  (E  +  iC0H)  =  Co(^  — v)(pc-j),  (29) 

and  separating  the  scalar,  vector,  and  bivector  parts, 
we  get 

0=|-Vq  (30a) 

□E  =  --VP-mo§,  (30b) 

e0  dt 

□H  =  V  x  j.  (30c) 

The  first  equation  is  the  continuity  condition[10], 
while  the  second  and  third  are  the  wave  equations 


3.3  Energy-Momentum 

The  product  of  the  electromagnetic  field  E  =  E  + 
i£oH  and  its  spatial  inverse  E$  =  — E  +  z£H  is 

E&  =  (E  +  *C0H)(-E  +  «CoH),  (33) 

which  may  be  expanded  as  [5,  6] 

EE*  =  — |E|2  -  Co|H|2  -  2CoE  x  H.  (34) 

This  expansion  lets  us  to  define  the  energy- 
momentum  cliffor  S/c: 


^=_IeoE#  =  t/+®, 

c  2  c 

where 

U=\(e0  |E|2+/r0|H|2), 

S^ExH 
c2  c2 


(35) 


(36a) 

(36b) 


are  the  energy  and  momentum  densities,  respectively. 
The  vector  S  =  E  x  H  is  the  Poynting  vector. 


4  Wave  Equation  in  Plasma 

4.1  Electron  Orbits 

If  we  assume  that  the  only  force  acting  on  the  charges 
is  significantly  due  to  the  electric  field  of  the  electro¬ 
magnetic  waves,  then  we  may  write  Newton’s  laws  of 
motion  as  [12] 

dv 

m7^2  =  «E-  (37) 


4 


16 


Note  that  we  assume  that  the  acceleration  of  the  elec¬ 
tron  is  mainly  due  to  the  electric  field  E  of  the  electro¬ 
magnetic  wave  and  contribution  of  the  wave’s  mag¬ 
netic  field  to  the  electron’s  motion  is  negligible. 

In  general,  the  electromagnetic  wave  is  ellipticaly 
polarized,  so  that  its  electric  field  may  be  expressed 
as:  [13] 

E  =  +  d+ij>)  (38) 

where 


a±  =  a±ele3a±, 


^±1  _  e±ie3(uit-kz) 


(39a) 

(39b) 


are  the  complex  amplitudes  and  the  wave  functions. 
Using  the  exponential  relation  in  Eq.  (13a),  Eq.  (38) 
may  be  expanded  as 


E  =  ei(a_  cos  (cot  —  kz  —  a_)  + 
a+  cos  (u>t  —  kz  +  a+))  + 
e2 (— a_  sin(wf  —  kz  —  a_)+ 

a+  sin(wf  —  kz  +  a+)).  (40) 

Though  Eq.  (40)  is  more  explicit  than  Eq.  (38),  the 
expression  involving  the  wave  functions  is  easier  to 
interpret  geometrically:  the  electric  field  E  is  a  sum 
of  two  counter-rotating  circular  motions.  The  cir¬ 
cles  have  magnitudes  a_  and  a+,  and  phases  a_  and 
a+.  For  fixed  position  z,  the  rotations  are  counter¬ 
clockwise  for  angular  frequency  u>  >  0  and  clockwise 
for  angular  frequency  ui  <  0.  For  fixed  time  t,  the 
rotations  are  clockwise  for  wave  number  k  >  0  and 
clockwise  for  wave  number  k  <  0. 

To  solve  equation  of  motion  in  Eq.  (37),  we  assume 
that  the  total  solution  r  may  be  written  as  solutions: 


r  =  rh  +  rp, 


(41) 


where  17  is  the  homogeneous  solution  and  rp  is  the 
particular  solution. 

The  homogeneous  solution  r h  satisfies  the  homo¬ 
geneous  equation 


m 


d2rh 
dt 2 


=  0. 


(42) 


From  this  we  see  that  the  solution  is  a  linear  motion: 


=  i-q  +  Vot,  (43) 


where  rg  and  Vg  are  unknown  vector  parameters  that 
still  have  to  be  determined  from  the  boundary  con¬ 
ditions. 

On  the  other  hand,  for  the  equation  for  the  partic¬ 
ular  solution  is 


d2rp 

dt,2 


qE. 


(44) 


If  we  assume  that  the  position  r  of  the  electron  is 
proportional  to  the  forcing  electric  held  E  in  Eq.  (38) , 
then 


rp  =  7E,  (45) 

where  7  is  a  constant.  Since  the  time  derivatives  of 
the  electric  held  E  in  Eq.  (38)  are 


<9E 

~dt 

c?2E 

dt2 


ei*e3w(— a,-^)  1  +  a+4>), 
ei(ie3w)2(d_'0'1  +  a+i >) 


(46a) 
-w2  E.  (46b) 


So  substituting  the  assumed  solution  in  Eq.  (45)  and 
using  the  expression  for  the  second  time  derivative  of 
the  electric  held  E  in  Eq.  (46b),  we  arrive  at 


g 

TOW2 


(47) 


Since  7  <  0,  then  the  particular  solution  rp  is  oppo¬ 
site  in  direction  to  the  forcing  electric  held  E. 

Substituting  the  results  for  17  and  rp  in  Eqs.  (43) 
and  (45)  in  Eq.  (41),  we  obtain 


V0t- 


rE, 


(48) 


Equation  (48)  is  the  total  solution  for  the  electron’s 
orbit  under  the  electric  held  E.  Notice  that  the  mo¬ 
tion  of  the  electron  is  a  sum  of  a  linear  motion  rg+Vgf 
and  an  elliptical  motion  opposite  to  the  electric  held 

E. 

Let  us  assume  that  at  t  =  0,  the  electron  is  sta¬ 
tionary,  so  that 


r  =  ro, 

(49a) 

v  =  0. 

(49b) 

With  these  boundary  conditions,  Eq.  (48) 
two  simultaneous  equations: 

leads  to 

ro  =  ro  ei(a_+a+), 

mujz 

(50a) 

0  =  Vq - — eiie3(— a_  +  a+). 

muj 

(50b) 
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Solving  for  ig  and  Vq,  we  get 

ro  =  r(H - ^ei(a_  +  a+),  (51a) 

mujz 

v'0  =  -^—e2(—d-  +  d+).  (51b) 

mu> 

Equations  (51a)  and  (51b)  are  the  expressions  for  the 
unknown  vector  parameters  of  the  charge’s  position 
in  time  in  terms  of  the  charge’s  initial  position  ro  and 
initial  velocity  Vo- 

Substituting  Eqs.  (51a)  and  (51b)  back  to  the  total 
solution  in  Ecp  (48),  we  obtain 

q  „  .  q  ,  „  „  . 

r  =  r0  H - -ei(a_  +  a+)  H - e2(-a_  +  a+)t 

mojz  muj 

- ^ei  (a-ip-1  +  a+ip).  (52) 

mujz 

Using  the  exponential  relation  in  Eq.  (13a),  Eq.  (52) 
expands  to 

r  =  ei  (xq  +  a_  cos  a_  +  a+  cos  a++ 

— ^-r-(((a_  sina_  —  a+  sina+)i 
mw2 

—  a_  cos  (tot  —  kz  —  a_) 

—  a+  cos  (u>t  —  kz  +  a+)))  + 
e2(yo  +  a-  sina_  +  a+  sina++ 

^  ((— a_  sina_  +  a+  sina+)t 

muiz 

+  a_  sin(a it  —  kz  —  a_) 

—  a+  sin(wt  —  kz  +  a+))).  (53) 

Equation  (53)  gives  the  orbit  of  the  electron  which  is 
initially  motionless  prior  to  its  interaction  with  an  el- 
liptically  polarized  electromagnetic  wave  of  frequency 
to. 

4.2  Wave  Equation 

4.2.1  Electric  Field 


Using  the  expression  for  the  total  solution  r  in 
Eq.  (52),  together  with  the  time  derivative  of  the 
electric  field  E  in  Eq.  (46a),  we  arrive  at 


j 


q  <9E 
raw2  dt 


(55) 


Notice  that  the  current  density  is  a  sum  of  an  con¬ 
stant  current  density  and  the  time  derivative  of  the 
electric  field  of  the  electromagnetic  wave. 

Substituting  the  expression  for  the  current  density 
j  in  Eq.  (55)  to  the  wave  equation  in  Eq.  (32),  we  get 


152e_v2f=  Nq2  d2g 

c2  dt2  ^  0  raw2  dt2 


(56) 


Using  the  expression  c  =  l/i//io£o,  we  may  rewrite 


Eq.  (56)  as 

,  v  '  ”  ”  ' 

n2  92E  2 

V2E  =  0 

c2  dt 2 

(57) 

where 

Nq2 

n=\l  2 

y  e0rmoz 

(58) 

is  the  plasma  refractive  index[l]. 

Equation  (58)  suggests  the  following  definition  for 
plasma  frequency: 


Nq 2 

Wo  —  - • 

eo  m 


(59) 


In  this  way,  the  refractive  index  n  in  Eq.  (58)  simpli¬ 
fies  to 

n = (6o) 

Notice  that  at  high  electromagnetic  wave  frequen¬ 
cies,  uj  »  wo,  the  ionosphere’s  refractive  index  ap¬ 
proaches  unity.  When  the  electromagnetic  frequen¬ 
cies  approaches  the  plasma  frequency,  w  — >  u>o,  the 
refractive  index  becomes  zero.  And  below  the  plasma 
frequency,  ui  <  wq,  the  refractive  index  is  imaginary. 


To  obtain  the  wave  equation  for  the  electric  field  E, 
we  assume  that  the  ionosphere  is  made  up  of  electrons 
of  charge  q.  If  N  is  the  number  of  charges  per  unit 
volume  and  these  charges  are  moving  with  velocity 
v,  then  the  current  density  j  is  given  by 


4.2.2  Magnetic  Field 

On  the  other  hand,  to  obtain  the  wave  equation  for 
the  magnetic  field  H,  we  take  the  spatial  derivative 
V  of  the  current  density  j  in  Eq.  (55)  to  obtain 


j 


(54) 


Vj 


(61) 


6 


18 


Separating  the  scalar  and  bivector  parts,  we  get 
N d 

v-j  = - %^V-E,  (62a) 

mcoz  ot 
N  c?  d 

Vxj  = - ^-VxE,  (62b) 

mo:z  dt 


equal  to  vacuum  permeability  no,  since  the  magnetic 
field  H  of  the  electromagnetic  wave  has  negligible  ef¬ 
fect  on  the  motion  of  the  charges  compared  to  that 
of  the  electric  field: 


A4  =  Mo- 


(67) 


after  factoring  out  i  in  the  second  equation. 

Using  Faraday’s  law  in  Eq.  (25c),  Eq.  (62b)  be¬ 
comes 


V  x  j  = 


Nq2n o  <92H 


(63) 


tow2  dt2 

Substituting  this  to  the  wave  equation  for  the  mag¬ 
netic  field  in  Eq.  (30c)  and  rearranging  the  terms,  we 
get 


n2  <92H 
c2  dt 2 


-  V2H  =  0, 


(64) 


where  we  used  the  definition  of  the  refractive  index 
n  in  Eq.  (58).  Equation  (64)  is  the  magnetic  field 
analog  of  the  plasma  wave  equation  for  the  electric 
field  E  in  Eq.  (57). 


For  the  permittivity  e,  we  know  that  the  wave  speed 
c/n  is  related  to  the  permittivity  and  permeability 
by 


c  1 

n  y/JIe' 


(68) 


Solving  for  the  permittivity  e  and  using  the  result  in 
Eq.  (67),  we  get 

e  =  n2e0,  (69) 


where  we  used  the  relation  c  =  l/-y/Moeo-  Thus,  the 
ratio  of  the  permittivity  of  the  plasma  and  vacuum 
is  the  square  of  the  refractive  index  n. 

Since  we  already  know  the  plasma  permeability  /i 
in  Eq.  (67)  and  the  plasma  permittivity  e  in  Eq.  (69), 
we  can  now  compute  the  plasma  impedance 


4.2.3  Electromagnetic  Field 

In  our  unperturbed  field  approximation,  we  used  the 
initial  electric  field  E  to  compute  the  current  density 
j  to  obtain  the  current  density  of  the  medium  as  a 
function  of  the  original  electric  field.  Clearly,  the 
original  electric  and  magnetic  fields  do  not  satisfy  the 
new  wave  equations  in  Eqs.  (57)  and  (64).  To  reflect 
this  change,  let  us  use  primed  variables  for  the  new 
electric  and  magnetic  fields  in  plasma,  so  that  the 
new  wave  equations  become: 


□'E'  =  0, 

(65a) 

□'H'  =  0, 

(65b) 

where 

□'  =  I^-V2 

c2  dt 2 

(66) 

is  the  d’Alembertian  operator  in  a 

medium  with  re- 

tractive  index  n. 

To  combine  the  electric  and  magnetic  wave  equa¬ 
tions  into  one,  we  need  to  compute  first  the  plasma 
impedance  Q  in  terms  of  the  plasma  electric  permit¬ 
tivity  e  and  plasma  magnetic  permeability  /i.  For 
the  magnetic  permeability  let  us  assume  that  it  is 


C  =  \/^=  “Co-  (70) 

V  e  n 

Thus,  the  ratio  of  the  impedances  of  the  plasma  and 
vacuum  is  equal  to  the  reciprocal  of  the  refractive 
index  n. 

We  may  now  construct  our  electromagnetic  field  E 
as 

E'  =  E'  +  *£H',  (71) 

so  that  the  two  wave  equations  in  Eq.  (65a)  and  (65b) 
combines  to 

D'CE'  +  i— H')  =  0.  (72) 

n 

Equation  (72)  is  now  a  homogeneous  electromagnetic 
wave  equation  in  a  plasma  with  refractive  index  n. 


4.3  EM  Wave 


But  there  are  difficulties.  If  we  use  the  expression  for 
the  current  density  j  in  Eq.  (55)  in  the  Ampere’s  law 
in  Eq.  (25b)  ,  we  get 


1  <9E' 
c  dt 


-CoVxH' 


-Co  Nq 


q  dE' 
tow2  dt 


,  (73) 
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where  we  already  replaced  the  fields  by  their  primed 
counterparts.  Rearranging  the  terms,  we  have 

^2  f)T?' 

^  -  CoV  x  H' = -CoA^v' .  (74) 

c  ot 

It  is  possible  that  this  extra  current  density  on  the 
right-hand  side  is  an  artifact  of  our  not  considering 
the  magnetic  field  of  the  electromagnetic  wave  in  the 
Lorentz  force  when  we  computed  the  orbit  of  the 
charges  Eq.  (37)  from  which  we  obtained  the  current 
density.  The  magnetic  field  would  provide  the  nec¬ 
essary  reaction  term  that  would  prevent  the  charges 
from  drifting  away.  In  this  paper,  we  shall  simply 
drop  the  Vq  term  in  the  Ampere’s  law: 

^-CoVxH'  =  °.  (75) 

c  Ot 

We  shall  also  assume  that  the  net  charge  density 
is  low.  That  is, 

p  «  0.  (76) 

This  may  be  a  good  approximation  in  the  ionosphere, 
since  the  electrons  are  knocked  out  from  the  atoms,  so 
that  the  effective  charge  density  is  zero.  The  positive 
ions  would  be  more  massive  than  the  electrons,  so 
that  the  positive  ions  would  move  much  more  slowly 
than  the  electrons.  Hence,  the  polarizability  of  the 
ionosphere  would  be  mainly  due  to  the  electrons. 

Equations  (74)  and  (76)  justify  the  following  ap¬ 
proximation: 


+  (E'  +  i^H')  =  °,  (77) 

\  c  ot  )  n 

which  is  the  homogeneous  Maxwell’s  equation  in 
plasma.  This  means  that  we  may  use  the  form 
of  the  electromagnetic  wave  in  vacuum  as  our  trial 
solution[13]: 

E'  =  E'  +  f-CoH'  =  e+(«-V’,~1  +  a'+V’O’  (78) 

n 


where 


e+  —  ei  +  i&2, 

—  g*e3  (ut—k'z) 


(79a) 

(79b) 


is  complex  basis  vector  and  the  plasma  wave  function. 

To  determine  the  wave  number  k'  in  Eq.  (79b), 
we  substitute  Eq.  (78)  back  to  the  wave  equation  in 
Eq.  (72),  we  obtain 


/  n2u 2 


E  =  0. 


(80) 


For  this  equation  to  be  satisfied,  the  first  factor  must 
be  zero,  so  that 

k'  =  nk ,  (81) 

where  k  =  uj/c.  Thus,  the  ratio  of  the  wave  number 
k'  in  plasma  and  that  of  k  in  vacuum  is  the  refractive 
index  n. 

Let  us  consider  three  cases:  n  >  0,  n  =  0,  and  n  is 
imaginary. 


4.3.1  Case  1:  n2  >  0 

If  n2  >  0,  then  the  electromagnetic  wave  E'  in 
Eq.  (78)  is  simply  an  clliptically  polarized  wave  which 
oscillates  in  both  time  and  space.  Separating  the  vec¬ 
tor  and  bivector  parts  of  Eq.  (78),  we  get 

E  =  +  a+x})'),  (82a) 

H  =  e2n(a'_xi)'~1  +  a+xpr).  (82b) 

These  may  be  expanded  as 

E'  =  ei(a'_  cos  (cot  —  nkz  —  a'_) 

+  a!+  cos  (ut  —  nkz  +  a+))  + 

G2(—d-  sin(wf  —  nkz  —  a'_) 

+  a!+  sin(wf  —  nkz  +  as.)) ,  (83a) 

H'  =  e2n(a'_  cos(wf  —  nkz  —  a'_) 

+  a'+  cos  (ut  —  nkz  +  oq_))  + 
e1n(o,_  sin(u;f  —  nkz  —  a!_) 

—  a!+  sin(wt  —  nkz  +  «+)) .  (83b) 

Notice  that  the  electric  and  magnetic  helds  are  both 
elliptically  polarized  since  they  are  both  a  sum  of  the 
same  counter-rotating.  The  only  difference  is  that 
the  two  helds  are  perpendicular  to  each  other. 
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4.3.2  Case  2:  n2  <  0 

And  if  n2  <0,  then  we  have  to  decide  on  which  unit 
imaginary  number  should  we  use,  because  there  are 
at  least  four  distinct  objects  which  square  to  —  1:  i, 
ie i,  ie 2,  and  *63.  From  the  form  of  the  wave  func¬ 
tion  ip+— '  in  Eq.  (79b),  we  see  that  the  proper  unit 
imaginary  number  should  be  ie 3,  so  that 

n  =  -b/e3  n-,.  (84) 

where 


is  the  magnitude  of  the  refractive  index. 

Using  the  expressions  in  Eqs.  (84)  and  (81),  the 
electromagnetic  wave  in  Eq.  (78)  becomes 

E'±  ±  ACoe3H'±  =  e+(a'_e~ie3Ute^kz 
\n\ 

+  a,+eie3Wte±Hfez),  (86) 

where  we  used  the  identities 

(ze3)_1  =  =  —  *e3,  (87a) 

te3 

(e3)2  =  1.  (87b) 

If  we  assume  that  H'±  is  in  the  xy— plane,  then 
e3  •  H±  =  0,  so  that 

e3H±  =  i(e3  x  Hj.)  =  — H-|-e3.  (88) 

is  a  bivector.  This  result  lets  us  separate  the  vector 
and  bivector  parts  of  Eq.  (86)  into 

E'±  =  e1(a,_e“ie3Wte=Fl"|fez  +  a'+eie3UJte±Hkz), 

(89a) 

C0H(j_  =  ±ei|  n\(a'_e~‘ie3Wte^n\kz  +  a'+eie3^e±Mkz), 

(89b) 

where  we  multiplied  the  bivector  part  by  ±e3  and 
used  the  identity 

e3*e2  —  e3e3el  =  ei-  (90) 


Expanding  the  exponentials,  Eqs.  (89a)  and  (89b) 
become 

E'±  =  ei (a'_e^n\kz  cos (ut.  —  c/_)+ 
a'+e±Mkz  cos(ujt  +  «'_))  + 
e2(-a'_e^nlkz  sin (wi  -  a'_)+ 

a'+e±|nlfez  sin (wt  +  a'_)),  (91a) 

CoH'±  =  ±ei\n\{a’_eT^kz  cos(wf  —  c/_)  + 
a'+e±|n|fcz  cos(wi  +  «'_))  + 

±  e2|?r|(— a'_eT^kz  sin(wt  —  a'_ ) 

a/+e±l"|fcz  sin  (wi  +  al)).  (91b) 

Notice  that  E±  ||  H±,  which  agrees  with  Budden. 

In  the  limit  as  z  »  l/|?r|fc,  the  terms  involving 
e~\n\kz  q;  so  that  Eqs.  (89a)  and  (89b)  reduces  to 

E'±  =  e1a(te|n|feze±ie3Wt,  (92a) 

CoH'±  =  ±e1|n|o^el"lfcze±ie3"t.  (92b) 

Expanding  the  exponentials,  we  get 

E'±  =  a±e^kz(e  1  cos  (u)t  ±  a'±)  ±  e2  sin  (u>t  ±  a±)), 

(93a) 

CoH'±  =  ±a±e'n\kz  (e  1  cos  (cot  ±  a'±)  ±  e2  sin  (cot  ±  o4))- 

(93b) 

In  other  words,  as  the  propagation  distance  z  be¬ 
comes  very  large,  the  electric  and  magnetic  fields 
become  circularly  polarized  and  their  amplitudes  in¬ 
crease  exponentially  with  distance. 

4.3.3  Case  3:  n2  =  0 

If  n2  — >  0+,  then  Eq.  (78)  reduces  to 

E'  =  E'  +  i 1CoH'  =  (ei  +  ie2){a!_e~ie3U>t  +  a',eie3Ut). 
n 

(94) 

Separating  the  scalar  and  bivector  parts, 

E'  =  e1(o/_e_*e3“t  +  a'+eie3Wt),  (95a) 
CoH'  =  0,  (95b) 
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so  that  only  the  electric  field  remains.  Expanding  the 
expression  for  the  electric  field,  we  arrive  at 

E'  =  ei(a_  cos {ut  —  a'_)  +  a+  cos (ut  +  $+))+ 
e-2(—a-  sin(wt  —  a_)  +  a+  sin(w<  +  4>+)). 

(96) 

Geometrically,  as  n2  —>  0+,  the  electric  field  E'  be¬ 
comes  a  sum  of  two  counter-rotating  circular  motions 
with  frequency  u>.  The  circles  have  radii  a_  and  a+, 
and  phase  angles  a'_  and  a'+ .  The  rotations  does  not 
change  phase  as  the  position  z  of  the  wave  varies. 

On  the  other  hand,  if  n2  — >  0“ ,  then  Eq.  (89a)  and 
(89b)  reduces  to 

E'±  =  e1(d/_e_ie3“t  +  a'+eie3Ut),  (97a) 
CoH'±  =  0,  (97b) 

which  are  the  same  equations  given  in  Eqs.  (95a)  and 
(95b). 

4.4  Wave  Energy-Momentum 


The  energy-momentum  of  an  electromagnetic  field  in 
plasma  is  is  analogous  to  that  in  vacuum  as  given  in 
Eq.  (35): 


-^eE'E^  =U'  +  ^/JTeS, 

(98) 

where 

U'  =  i(eE'2+AcH'2), 

(99a) 

lieS'  =  /jeE'  x  H', 

(99b) 

In  terms  of  the  refractive  index  n  and  the  speed  of 
light  in  vacuum  c,  Eq.  (98)  may  be  written  as 

-\eE'E^  =  U'+  -S', 

2  c 

(100) 

where 

U'  =  \{e0n2  E'2  +  Ai0H'2), 

(101a) 

^S'  =  ^E'  x  H',  (101b) 

where  we  used  the  relations  /i  =  /iy  and  e  =  rc2e o 
in  Eqs.  (67)  and  (69),  respectively.  Eqs.  (101a)  and 


(101b)  are  the  energy  and  momentum  densities  of  the 
electromagnetic  field  in  plasma. 

Let  us  compute  the  energy-momentum  of  the  elec¬ 
tromagnetic  wave  E'  in  Eq.  (78),  by  considering  three 
cases:  n2  >  0,  n2  =  0,  and  n2  <  0. 

4.4.1  Case  1:  n2  >  0 

For  n2  >  0,  the  energy  momentum  of  the  electromag¬ 
netic  wave  E'  in  Eq.  (78)  is 

U'  +  —S'  =  —  ^eo«2e+(d'_^/_1  + 
c  2 

e^(a,_'i/),_1  +  (102) 

since 

a],  =  a±,  (103a) 

(^'±1)t=^'±ij  (103b) 

because  scalar-bivector  cliffors  are  invariant  under 
spatial  inversion. 

From  the  relations  in  Eqs.  (13a)  and  (13b),  we  have 

e+i/^1  =  '0=Fle_|_ ,  (104a) 

e+a±  =  a±e+.  (104b) 

Hence,  Eq.  (102)  becomes 

U'  +  -S'  =  -;^-eo?^2e+e1|_(d'_!V’,  + 
c  2 

(a,_^,_1  +  d'+V>'),  (105) 

Distributing  the  terms  in  Eq.  (105)  and  using  the 
identity 

e+e(|_  =  —2(1  +  e3)  (106) 

we  obtain 

U'  +  -S'  =  )-e0n2(l  +  e3)(o2  +  ci2+ 
c  2 

+  ala'+xi)'2  +  d/_Ja,_-i/),_1).  (107) 

Using  the  definitions  of  Eq.  (39a)  and  (79b), 
Eq.  (107)  becomes 

U'  +  —S'  =  (1  +  e3)eo  n2(a12  +  a1?  +  2  a'_a'+x 
c 

cos(2 (uit  —  nkz)  +  a'+ —  (108) 
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Separating  the  scalar  and  vector  parts  of  Eq.  (108), 
we  get 

U'  =  eo  n2{a'2  +  a'2  +  2  a!_a'+  x 

cos(2  {u>t  —  nkz)  +  a'+ —  a'_)),  (109a) 


cos(2  (cut  —  nkz)  +  a'+ —  a'_)).  (109b) 

That  is,  both  energy  density  U'  and  the  momentum 
density  n2S'/c2  are  fluctuating  in  space  and  time  at 
twice  the  frequency  of  the  incident  electromagnetic 
wave.  Note  that  the  Poyning  vector  S'  is  proportional 
to  the  refractive  index  n: 

S'  =  e0n(a12  +  a12  +  2  a'_a'+  x 

cos(2(w<  —  nkz)  +  a'+  —  ol_)).  (HO) 

Averaging  Eq.  (108)  over  the  optical  period  r  = 
27 r/w, 

(U1  H — S')T  =  (1  +  e3)e3n2  (a'z  +  o,'2),  (HI) 
c 

and  separating  the  energy  and  momentum  compo¬ 
nents,  we  arrive  at 

( U')T  =  eon2  (a12  +  a'2),  (112a) 

77^  77^ 

^<S')T  =  e3eo-(a'_2  +  a'2).  (112b) 

c  c 

Thus,  the  average  energy  is  proportional  to  the 
square  of  the  refractive  index  n  while  the  average 
momentum  is  proportional  to  the  cube.  In  partic¬ 
ular,  the  average  Poynting  vector  is  proportional  to 
the  refractive  index  n: 

(S')T  =  e3e0nc{a12  +  a12).  (113) 

4.4.2  Case  2:  n2  <  0 

If  n2  <  0,  we  use  the  electromagnetic  wave  expression 
in  Eq.  (86).  Its  corresponding  energy- momentum 
equation  is 

^  _2  £,/  jVf  ^  (  X'  —ie^Lot  ^flnlkz 

--e0n  E±EJ_  = --e0n  e+{a_e  3  e+l  1 

+  a'+eie3Ute±Mkz) 
e\_(a'_e~ie3Ute^kz 
+  d'+eie3Ute±Mkz).  (114) 


Using  the  identities  in  Eqs.  (104a),  (104b),  and  (106), 
Eq.  (114)  becomes 

-\eon2E’±E'l  =  (1  +  e3)e0n2{a'*eie3Ute^kz 
+  d'*e-ie3Ute±ln\kz) 

(a'  e~ie3Ut e^n^kz 
+  a'+e~ie30Jte±\n\kz).  (115) 

Distributing  the  terms  in  Eq.  (115)  yields 

-  *e0  n2E'±E’l  =  (1  +  e3)e0n2  (a'2e*2Wkz  +  a,2e±2^kz 

a'*a'_e~2ie3Ut  +  d'*a+e2ie3Ut). 

(116) 

Equation  (116)  is  the  expansion  of  the  right  side  of 
Eq.  (114). 

To  expand  the  left  side  of  Eq.  (114),  we  use  the 
definition  of  the  electromagnetic  field  E '  for  n 2  <  0 
in  Eq.  (86),  together  with  the  : 

-\e0n2E'±E'l  =  -±en2  (V±  ±  ^Coe3H'±) 

(-E'±  ±  |^Coe3H'±)  ,  (117) 

since 

(e3H'±)t  =  ejH'2  =  e3H'±  (118) 

and  the  spatial  inverse  of  a  vector  is  its  negative. 
Distributing  the  terms  in  Eq.  (117), 

eon2E'±E'l  =  -\e0n2{- E'|  +  ^e3H'±e3H'± 

±^(E'±e3H'±^e3H'±E'±)), 

(119) 

and  using  the  perpendicularity  assumption  in 
Eq.  (88),  we  get 

^  2  Pi/  zVf  _  1  t12('TT'2  i  Co  tt/2 

~2e°n  L±h±  -  2  0  ^E±  +  |u)2H± 

i  ■j~j'e3(E,±H(t  +  H(j_E(j_)).  (120) 
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From  Eqs.  (89a)  and  (89b),  we  see  that 

E'+=±— H±.  (121) 

n 

Using  this  equivalence  in  Eq.  (120),  we  obtain 

-^e0n2E'±E'l  =  e0n2(l  +  e3)E'|,  (122) 

which  is  a  simple  expression  in  terms  of  E(jr.  Com¬ 
paring  Eq.  (122)  with  Eqs.  (99a)  and  (99b)  leads  to 
the  relation 

~e0n2E'±E'l  =  U±  +  s'±,  (123) 

where 

4  =  e0n2e3E'±2  (124) 

is  proportional  to  momentum  density.  Thus,  even  if 
the  wave’s  Poynting  vector  is  zero, 

S'±  =  E±  x  H'±  =  0,  (125) 

because  E'±  and  H'±  are  parallel  by  Eq.  (121)  as  noted 
by  Budclen,  the  electromagnetic  wave  still  possess  a 
quantity  s'±  with  units  of  energy,  which  we  interpret 
as  proportional  to  momentum  density  .  This  is  a  new 
result. 

Equating  Eqs.  (123)  with  (116),  we  get 

U'±  +  s'±  =  (1  +  e3)e0n2(a,^e=F2|n|fez  +  a'2e±2^kz 
a!*a!_e-ie32ut  +  a'la+eie32ut) .  (126) 

Expanding  the  exponentials, 

TTf  i  _/  / -i  |  \  2  /  /2  ^E2\n\kz  ,  /2  „±.2\n\kz 

U±  +  s±  =  (1  +  e3)e0n  (a_eT  1  1  +  a+e  1  1 

+  2a'+a'_  cos(2w<  +  a'+  —  a'_)),  (127) 

and  separating  the  scalar  and  vector  parts,  we  arrive 
at 

U'±  =  e0n2{a'2eT2Mkz  +  a,2e±2^kz 

2a+a-  cos(2wf  +  a'+  —  a_),  (128a) 

_  _  ,  2  /  /2  „^f2\n\kz  ,  /2  „±2\n\kz 

s±  =  e3eon  (a_e^  1  1  +  a+e  1  1 

2a+a_  cos(2wf  +  a'+  —  a_).  (128b) 

Equations  (128a)  and  (128b)  are  the  desired  rela¬ 
tions  for  energy  and  momentum  densities  for  the  case 


n2  <  0.  ince  n2  <  0,  then  both  the  energy  density 
U±  would  be  negative,  and  s±  (a  quantity  we  claim 
to  be  proportional  to  momentum  density)  would  be 
pointing  in  the  direction  opposite  to  e3,  the  direction 
of  motion  and  momentum  of  the  wave  if  n2  >  0. 

Taking  the  time  average  of  the  energy  and  momen¬ 
tum  relations  in  Eqs.  (128a)  and  (128b)  over  the  pe¬ 
riod  t  —  2tt/lo,  we  obtain 

U'±  =  e0nV2eT2|”|fez  +  a,2e±2^kz),  (129a) 

s'±  =  e0n2(a,2eT2|n|fcz  +  a,2e±2^kz).  (129b) 

Depending  on  the  choice  of  the  sign  of  the  refractive 
index  n  =  ±*e3|n|,  either  a'_  or  a'+  will  dominate  as 
z  »  l/\n\k. 

4.4.3  Case  3:  n 2  =  0 

For  n2  0+,  the  electromagnetic  wave  E'  is  given 
in  Ecj.  (94).  To  find  its  energy- momentum,  we  set 
n  =  0  in  Eq.  (108)  to  obtain 

2 
77  ^ 

U'  =  — S'  =  0.  (130) 

c2 

On  the  other  hand  for  n2  — >  0_,  the  electromag¬ 
netic  wave  E±  is  given  in  (86).  To  find  its  energy- 
momentum,  we  set  n  =  0  in  Eq.  (126)  to  get 

U±  =  s±  =  0.  (131) 

Regardless  of  the  direction  of  evaluation  of  the  limit 
of  \n\  —>  0,  we  get  the  same  conclusion:  the  electro¬ 
magnetic  wave  possesses  neither  energy  nor  momen¬ 
tum. 

5  Conclusions 

In  this  paper  we  used  geometric  algebra  for  comput¬ 
ing  the  energy  and  momentum  densities  of  an  cllip- 
tically  polarized  wave  in  plasma  for  different  cases  of 
refractive  indices:  n2  >  0,  n2  <  0,  and  n 2  =  0. 

In  Section  2,  we  began  by  introducing  Clifford  (ge¬ 
ometric)  algebra  CZ3jo  and  showed  that  the  direct 
product  of  two  vectors  can  be  expressed  as  a  sum  of 
their  dot  product  and  their  imaginary  cross  product, 
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with  the  unit  imaginary  number  as  the  unit  trivec¬ 
tor,  which  is  a  product  of  three  basis  vectors  along 
the  x— ,  y—.  and  z— directions.  We  showed  that  the 
exponential  of  an  imaginary  vector  may  be  used  as 
a  vector  rotation  operator,  for  vector  operands  that 
are  perpendicular  to  the  vector  argument  of  the  ex¬ 
ponential. 

In  Section  3,  we  summarized  the  fundamental 
equations  of  Electrodynamics.  The  first  equation  is 
the  Maxwell’s  equation,  which  states  that  the  space- 
time  derivative  of  the  electromagnetic  field  is  pro¬ 
portional  to  the  charge-current  density.  The  second 
equation  is  the  wave  equation,  which  is  obtained  by 
taking  the  derivative  of  the  Maxwell’s  equation  with 
respect  to  the  spatial  inverse  of  the  spacetime  deriva¬ 
tive  operator.  Here  the  current  density  is  a  sum  of  the 
electric  field  vector  and  the  imaginary  magnetic  field 
vector.  We  showed  that  the  product  of  the  electro¬ 
magnetic  field  with  its  spatial  inverse  is  proportional 
to  the  energy-momentum  of  the  electromagnetic  field. 

In  Section  4,  we  deduced  the  orbital  motion  of  a 
charge  subject  to  the  electric  field  of  a  circularly  po¬ 
larized  electromagnetic  wave.  We  showed  that  the 
motion  of  the  charge  is  a  superposition  of  a  linear 
motion  and  an  elliptical  motion,  with  the  latter  op¬ 
posite  to  the  forcing  electric  field.  We  took  the  time 
derivative  of  the  charge’s  position  vector  to  obtain  the 
current  density.  We  showed  that  the  current  density 
has  two  components:  a  constant  current  density  and 
a  time-varying  current  density  proportional  to  the 
time-derivative  of  the  electric  field  of  the  electromag¬ 
netic  wave.  We  dropped  the  constant  current  density 
term,  because  we  felt  that  this  is  only  an  artifact  of 
our  not  also  considering  the  magnetic  field  of  the  elec¬ 
tromagnetic  wave  in  the  Lorentz  force  on  the  charge. 
We  also  assumed  that  the  charge  density  is  nearly 
zero.  In  this  way,  we  showed  that  we  can  define  new 
electric  and  magnetic  fields  in  plasma  which  obeys 
a  homogeneous  wave  equation  and  a  homogeneous 
Maxwell’s  equation,  with  refractive  indices  computed 
from  the  proportionality  of  the  time  derivative  of  the 
current  density  to  the  second  time  derivative  of  the 
forcing  electric  field  of  the  electromagnetic  wave. 

We  assumed  an  electromagnetic  wave  similar  in 
form  to  that  in  vacuum  for  different  cases  of  the  re¬ 
fractive  index  n.  For  n2  >  0,  we  showed  that  the 


electric  and  magnetic  fields  are  perpendicular;  while 
for  n2  <  0,  the  two  fields  are  parallel.  For  n2  >  0,  the 
energy-momentum  of  the  wave  is  well-defined,  which 
is  a  sum  of  the  energy  density  and  the  Poynting  vec¬ 
tor.  For  n2  <  0,  the  energy  density  is  easy  to  inter¬ 
pret,  since  it  is  just  proportional  to  the  square  of  the 
wave’s  electric  field.  The  momentum  density  is  more 
difficult  to  interpret  since  the  Poynting  vector  is  zero, 
because  the  electric  and  magnetic  fields  are  parallel, 
as  noted  by  Bud  den.  Nevertheless,  we  showed  that 
we  can  construct  a  new  quantity  with  units  of  energy 
that  points  along  the  direction  of  propagation  or  rota¬ 
tion  of  the  wave;  we  claim  that  this  is  proportional  to 
momentum  density.  But  since  the  energy  density  and 
the  new  momentum  density  are  both  proportional  to 
n2  <  0,  then  the  energy  density  is  negative  and  the 
new  momentum  density  is  opposite  to  the  direction 
of  propagation  of  the  wave  for  the  case  n2  >  0. 

There  are  many  results  in  the  paper  which  may  not 
be  physically  correct,  such  as  drifting  charges  under 
polarized  light,  negative  energy  densities  of  electro¬ 
magnetic  waves,  and  the  exponential  increase  of  wave 
amplitude  as  function  of  propagation  distance.  We 
hope  to  remove  all  these  problems  in  our  succeeding 
paper  where  we  compute  the  refractive  index  from 
the  the  orbit  of  charges  under  the  action  of  both  the 
electric  and  magnetic  field  components  of  the  ellipti- 
cally  polarized  electromagnetic  waves.  We  shall  also 
try  to  extend  our  geometric  algebra  formalism  to  the 
reformulation  of  Appleton’s  magnetoionic  theory  of 
the  ionosphere  and  the  meaning  of  the  complex  re¬ 
fractive  index. 
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